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I. INTRODUCTION 

The relevance of a coherent relativistic description 
of few-body systems is now well recognized in parti- 
cle as well as in nuclear physics. Concerning particle 
physics, a relativistic formalism is necessary for the 
understanding of the various components of the nu- 
cleon or pion state vectors in terms of valence quarks, 
giuons, and sea quarks, as revealed, for instance, in 
exclusive reactions at very high momentum transfer. 
The need for a coherent relativistic approach to few- 
body systems has also become clear in nuclear physics 
in order to check the validity of the standard descrip- 
tion of the microscopic structure of nuclei in terms 
of correlated pion exchanges between nucleons within 
the general framework of chiral perturbation theory. 
In this case, electromagnetic interactions play a cen- 
tral role in "seeing" pion exchanges in nuclei. 

In the non-relativistic limit (when the speed of light 
c goes to infinity) a system of particles is described by 
its wave function defined at fixed moments of time or, 
in other words, on the plane t ~ const, and its time 
evolution is governed by the Schrodinger equation, 
once the Hamiltonian of the system is known. Rel- 
ativistic description admits some freedom in choosing 
the space-like hyper-surface on which the state vector 
is defined A possible choice is to take, for this 
purpose, the same plane t = const (the so-called "in- 
stant" form of dynamics). This is however not very 
well suited for relativistic systems, since this plane is 



not invariant under Lorentz boosts. It is much more 
preferable to use Light-Front Dynamics (LFD) which 
is of particular interest among various approaches ap- 
plied so far to study relativistic systems. In the stan- 
dard version of LFD, the state vector is defined on the 
plane ^ + ■§ = [Ij] , invariant with respect to Lorentz 
boosts along the z axis. 

Advantages of using LFD to describe physical sys- 
tems are well known. The main one concerns the 
structure of the vacuum. Because of kinematical con- 
straints, the plus-component p+ = p° -f (we take 
hereafter c = 1) of the four- momentum p of any par- 
ticle state, both real and virtual, is always positive 
or null. This implies that the vacuum state coincides 
with the free vacuum, and all intermediate states re- 
sult from fluctuations of the physical system. One 
can thus construct any physical system in terms of 
combinations of free fields, i.e. the state vector is de- 
composed in a series of Fock sectors with an increasing 
number of constituents. This enables a systematic cal- 
culation of state vectors of physical systems and their 
observables. 

Note that the triviality of the vacuum in LFD, men- 
tioned above, docs not prevent from non-perturbative 
zero-mode contributions (states with p+ = 0, some- 
times called the "vacuum sector") to field operators, 
when physical systems with spontaneous symmetry 
breaking are considered 0]. An application to the 0* 
model in 1 -|- 1 dimension has been done in Ref. 0]. 

While the Fock decomposition is non-perturbative, 
it is only meaningful if it converges rapidly. One way 
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to look at this convergence for a simple but never- 
theless physically relevant system is to investigate, 
within LFD, the Wick-Cutkosky model: a system of 
two scalar particles of mass m interacting by the ex- 
change of a massless scalar particle. Independently, 
the same system can be considered within the four- 
dimensional Fcynman approach by solving the Bcthc- 
Salpeter equation in the ladder approximation which 
includes exchanges of an infinite number of scalar 
bosons in the intermediate state. Comparing the re- 
sults of both calculations [1] , we can see that the two- 
and three-body components of the state vector repre- 
sent as much as 90% of its norm, for m = 1 GeV and a 
coupling constant of 2t: which gives the maximal bind- 
ing. Such a simple test shows that even in the worst 
case (a large coupling constant and the exchange of a 
boson of zero mass) the Fock decomposition is mean- 
ingful and may converge rapidly. This however should 
be analyzed in more realistic calculations. 

The decomposition of the state vector of any physi- 
cal system in terms of Fock sectors on the Light Front 
(LF) enables a very intuitive interpretation of the 
physical state, since each Fock sector is reminiscent 
of a non-relativistic many-body wave function. 

The standard version of LFD has however a serious 
drawback, since the equation of the LF plane t + z = 
is not invariant under spatial rotations. As we shall 
see later on, the breaking of the rotational invariance 
has many important consequences as far as the con- 
struction of bound states with definite angular mo- 
mentum is concerned, or in the calculation of electro- 
magnetic amplitudes. 

To avoid such an unpleasant feature of standard 
LFD, we shall use below the Covariant formulation of 
LFD (CLFD) 0,11], which provides a simple, practi- 
cal, and very powerful tool in order to describe phys- 
ical systems as well as their electromagnetic ampli- 
tudes. In this formulation, the state vector is de- 
fined on the plane characterized by the invariant equa- 
tion Lo-x = 0, where uj is an arbitrary light-like four- 
vector with = 0. The standard LFD on the plane 
t + z ~ is recovered by considering the particular 
choice lli = (1,0,0,-1). The covariance of our ap- 
proach is caused by the invariance of the LF plane 
equation uj-x ~ under any Lorentz transformation 
of both Lo and x. This implies in particular that lo 
cannot be kept the same in any reference frame, as it 
takes place in the standard formulation of LFD with 
u;= (1,0,0,-1). 

There is of course equivalence, in principle, between 
the standard and covariant forms of LFD. Within the 
same approximation (or for exact calculations) CLFD 



reproduces the results of standard LFD as a partic- 
ular case. The physical observables should coincide 
in both approaches, though their derivation in CLFD 
in most cases is much simpler and more transparent. 
The relation between CLFD and standard LFD re- 
minds that between the Feynman graph technique and 
old-fashioned perturbation theory. 

CLFD has first been used to investigate the gen- 
eral structure of few-body systems and their electro- 
magnetic observables in the tree approximation (sec 
Ref. @ for a review). If one wants to go beyond this 
phenomenological analysis, one has to be able to cal- 
culate the state vector of a physical system from a 
given Hamiltonian in a non-perturbative framework. 

Consider, as an example, a system composed of in- 
teracting fermion and bosons. In the simple two-body 
Fock space truncation, the physical fermion state vec- 
tor is represented as a sum of two sectors: the one sin- 
gle fermion state and the one fermion plus one boson 
state. The fermion propagator is thus given, in the 
chain approximation, by the contributions indicated 
in Fig. [ija). It is non-perturbative in the sense that 
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FIG. 1: Fermion propagator in the chain approximation, 
within the two-body Fock space truncation (a) and an ir- 
reducible contribution to the perturbative self-energy (b). 
Solid and wavy lines correspond to fermions and bosons, 
respectively. 

it involves contributions to all orders in the coupling 
constant g, but approximate, since it incorporates at 
most two particles in the intermediate states. It is well 
known that this infinite series can be summed up in 
terms of the (perturbative) self-energy S(p) of order 
g^, as indicated in Fig.[T]Jb). In this two-body trunca- 
tion, the equivalence between the LF fermion propa- 
gator (calculated in CLFD) and the two-point Green's 
function (calculated in the Feynman four-dimensional 
approach) has been shown to occur very naturally to 
all orders in 5 Q. 

The fermion propagator enters into the expression 
for the observable fermion-boson scattering ampli- 
tude. This amplitude must have a pole, in the s- 
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channel; at s = m^. To ensure such a property, a Mass 
Counterterm (MC) must be added to the self-energy. 
Besides that, the coupling constant coming into the 
vertices of the diagrams in Fig. [T] can not be identi- 
fied a priori with the physically observed quantity, but 
should be treated as some bare (non-renormalized) pa- 
rameter. In order to calculate physical observables, 
the boson-fermion Bare Coupling Constant (BCC) as 
well as the MC should be expressed in terms of the 
physical coupling constant and the particle masses. 
This has been done, for the two-body Fock state trun- 
cation in CLFD, in Ref. 0. However, a general renor- 
malization scheme needed to determine the MC and 
the BCC for the most general case of Fock space trun- 
cation has not been proposed yet. 

Already at the level of the two-body Fock space 
truncation, one has to deal with loop diagrams [like 
the self-energy contribution shown in Fig. [1] (b)]. 
Their amplitudes diverge for high internal momenta. 
The implementation of any rcnormalization scheme 
essentially depends on the way of regularization of 
divergent amplitudes. This is indeed a non-trivial 
task, as it has been already mentioned in various con- 
texts d, Q . The regularization of amplitudes in LFD 
by traditional cutoffs imposed on the transverse and 
longitudinal components of particle momenta, for in- 
stance, corresponds to restricting the integration vol- 
ume by a rotationally non-invariant domain. The reg- 
ularized amplitudes depend therefore not only on the 
size of this domain (i.e., on the cutoff values), but also 
on its orientation determined by the orientation of the 
LF plane. 

Another source of violation of rotational invariance 
is the Fock space truncation itself. As a consequence, 
the number and the structure of the counterterms 
needed to renormalize the theory depend on the LF 
plane orientation as well. CLFD allows us to param- 
eterize the latter dependence in a very transparent 
form, through the four-vector lo. Moreover, the co- 
variant formulation of the approach is mandatory in 
order to define what arc the physical parameters of 
the theory (and hence to be able to renormalize the 
latter), since it enables an explicit separation of any 
spurious contributions depending on lo. This is the 
case, for instance, for the two-body wave function, as 
we shall see in Sec. |lll 

Following the analysis of Ref. Q, we choose the 
Pauli-Villars (PV) regularization scheme in order to 
impart mathematical sense to divergent amplitudes. 
This scheme also preserves rotational invariance, as 
well as other important symmetries like gauge invari- 
ance. Though the PV regularization was developed 



initially for the four-dimensional Fcynman approach, 
it can be easily implemented into the LFD calculating 
machinery by simply introducing additional fictitious 
PV fields [ll. 

The rcnormalization procedure must ensure that 
physical results do not depend on the regularization 
parameters. Besides that, it should be, first, non- 
perturbative and, second, consistent with the trunca- 
tion of the Fock decomposition in the sense that it 
should not leave any divergences uncancelled. 

Let us look, for example, at the rcnormalization of 
the fermion propagator in the second order of pertur- 
bation theory. There exist three contributions to the 
physical fermion propagator, as indicated in Fig. [21 
These are, from left to right, the free propagator, 
the self-energy contribution S(p), and the contribu- 
tion from the MC 5m. The sum of these three items 
should be equal, at = m^, to the free propagator. 
This fixes 5m = — E(p) at ^ = m. 




FIG. 2: Renormalization of the fermion propagator in the 
second order of perturbation theory. 

As we can see from Fig. [2l Fock sectors with differ- 
ent number of constituents are intimately connected 
through the renormalization condition: the contribu- 
tion of the MC (the last diagram in Fig. [5]) corre- 
sponds to the one-body Fock sector (a single fermion). 
It should however be opposite, at ^ = to, to the two- 
body (one fermion plus one boson) Fock sector contri- 
bution given by the second diagram in Fig. [51 in order 
to cancel its divergence. This means that any MC or, 
more generally, any bare parameter, should be asso- 
ciated with the number of particles in a given Fock 
sector. In other words, all MC's and bare parameters 
must depend on the Fock sector under consideration. 
This is a necessary condition. 

Several attempts have already been made to ad- 
dress the problem of non-perturbative renormaliza- 
tion in the standard formulation of LFD, either in the 
Yukawa model (a fermion coupled to scalar bosons) or 
in QED, using various regularization schemes. Early 
calculations were performed with a momentum cut- 
off for the Yukawa model ^ and for QED As 
shown in Ref. [1] , the use of such a cut-off implies to 
consider specific counterterms depending on the LF 
plane orientation. Moreover, the absence of Fock sec- 
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tor dependent counterterms and BCC's prevents any 
calculation to converge properly. 

The use of PV fields to regulate the amplitudes 
has first been advocated in Refs. [l^l (with three 
PV bosons) and [l^ (with three PV fermions), for 
the Yukawa model and QED, respectively. These 
calculations suffer however from the lack of a non- 
perturbative procedure to determine the parameters 
of the PV fields, as well as from an incorrect chiral 
limit. Again, no Fock sector dependent counterterms 
were considered, which left divergences uncancelled. 
In particular, this prevents the two-body calculation 
of QED to reproduce the well known radiative cor- 
rection to the anomalous magnetic moment of the 
electron (the Schwinger correction). We shall see in 
Sec. llVl how it arises naturally in our scheme. 

Most recent calculations in the Yukawa model with 
the two- [ij] and three-body Fock space trun- 
cations used simultaneously a PV fermion and a PV 
boson to regulate the amplitudes. This regularization 
procedure is adequate to preserve rotational invari- 
ance, at least for the two-body truncation, according 
to the analysis of Ref. Q . However these two calcula- 
tions are plagued with uncancelled divergences. 

The dependence of the counterterms on the Fock 
sectors has been first suggested in Ref. [l6j in the con- 
text of a simple calculation within the two-body Fock 
space truncation. This idea has however never been 
formulated as a coherent renormalization scheme. 

The main aim of the present article is to develop 
such renormalization scheme. We propose a complete 
and coherent strategy to determine the counterterms 
and the bare parameters in LFD calculations with a 
Fock space truncation to any order. A prelimi nary 
account of such a scheme was presented in Ref. ■ 
We conjecture that this renormalization scheme is also 
sufficient to avoid any uncancelled divergences in any 
order of the Fock space truncation, provided appro- 
priate counterterms necessary to recover rotational in- 
variance (if needed) are taken into account. We shall 
demonstrate below that this is indeed the case for 
some model and realistic physical systems, within the 
two- and three-body Fock space truncations. 

The plan of our paper is the following. In Sec. |lTl 
we recall the main features of the description of bound 
state systems in CLFD, taking the Yukawa model as 
an example. We expose in Sec. lIIII our new systematic 
renormalization scheme in CLFD calculations with 
Fock space truncation. Applications of this scheme to 
particular physical systems — to the Yukawa model 
and QED — within the two-body Fock space trunca- 
tion (Sec. HV]) , and to a purely scalar model for the 



three-body truncation (Sec. |V|) are then considered. 
We present our concluding remarks and outline possi- 
ble perspectives in Sec. lVIl Some technical derivations 
are given in Appendices. 

II. DESCRIPTION OF PHYSICAL SYSTEMS 
IN THE COVARIANT FORMULATION OF 
LIGHT-FRONT DYNAMICS 

In order to show how our renormalization scheme 
should be applied to the analysis of physical systems, 
we shall consider in the following study the Yukawa 
model, i.e. a physical fermion composed of a bare 
fermion coupled to scalar bosons. This system is inter- 
esting from several points of view. It is not as simple 
as a super-renormalizable purely scalar model, while 
it has many similarities with QED in the Feynman 
gauge, at least for the case of the simple two-body 
truncation. It is thus easy to extend our results, as 
shown in Sec. IIVI 

A. The Yukawa model. Construction of the 
hght-front interaction Hamiltonian 

The Lagrangian describing a system of interacting 
spin-1/2 fermion and scalar boson fields, taking into 
account the mass renormalization of the fermion, is 

C^Cf + Cb+Cfb. (1) 

where the three terms on the r.-h.s. are, respectively, 
the fermion, boson, and interaction parts, 

Cf = i^-i^dy-^ -m^-^, (2a) 

Cb = i [a,$9''<i> - , (2b) 

-Cfb = .90^** + <5mi'^'. (2c) 

Here ^E* = '^{x) and $ = $(a;) are the Heisenberg 
fermion and boson field operators, 170 is the BCC, anal- 
ogous to the bare charge eo in QED, m is the physical 
fermion mass, /i is the physical boson mass, and 5m 
is the fermion MC. 

As already advocated in Ref. 0], it is more appro- 
priate and physically sounded to construct Fock sec- 
tors in terms of free fields corresponding to particles 
with their physical masses. In that case, one does 
not have to determine the fermion bare mass mo but 
rather a MC 5m = mo — m, as it is well known p^ . 
This choice of the renormalization procedure for the 
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fcrmion mass is the only way to keep the basis con- 
structed from free fields to be the same in all Fock 
sectors. In our renormalization scheme the bare pa- 
rameters like mo depend on the Fock sector in which 
they appear. If one assigned the bare mass mo to the 
free fermion field, the latter would be different in dif- 
ferent Fock sectors. Taking the free fermion field with 
the physical mass m, on the contrary, fixes it once and 
for all, while dependence of renormalization parame- 
ters on the Fock sectors is carried over to the MC Sm. 
Moreover, one may hope that the Fock state expan- 
sion may converge more rapidly with the use of a fixed 
physical mass as compared to a (divergent) bare mass. 

Working in LFD, we have to deal with Hamiltoni- 
ans, rather than Lagrangians. Moreover, since wc use 
Fock expansions in terms of free fields, the Hamil- 
tonian must be also expressed through them (i. e. 
taken in Schrodingcr or interaction representation). 
The general procedure of deriving CLFD Hamiltoni- 
ans from Lagrangians is exposed in Ref. 0. First, one 
should construct the energy-momentum tensor 

where Yi denotes either 5* or '5, or $, the sum run- 
ning over all the fields, and the LF four-momentum 
operator 

Pp = \j da''{x)Q,„ (4) 

where the integration is performed on the three- 
dimensional space element orthogonal to the "time" 
direction (the role of time is played in CLFD by the in- 
variant combination oj-x). The four-momentum oper- 
ator should then be expressed through the free fields, 
taking into account constraints imposed on the field 
components by the equations of motion. The corre- 
sponding operator Pp can be represented as the sum 

Pp = p^") + p;"*, (5) 

where the two terms on the r.-h.s. are, respectively, 
the free (i.e. independent of the coupling constant 
and countcrterms) and interaction parts of the four- 
momentum. The operator P™* is related to the inter- 
action Hamiltonian iJ™* (x) by 

P;"* ^Wp j W\x) 5{lo-x) d'^x. (6) 

The calculations performed in Ref. [7] showed that 
the interaction Hamiltonian for the Yukawa model in- 
cludes also a set of so-called contact (or instantaneous) 



terms which explicitly depend on the LF plane orien- 
tation and essentially complicate calculations, both 
pcrturbative and non-perturbative. 

Wc shall use hereafter the PV rcgularization which 
not only maintains rotational invariance, but also kills 
the contact terms, as will be demonstrated below. The 
PV scheme can be easily implemented into the La- 
grangian [l^ by introducing additional fields (we will 
call them PV fields or PV particles), having negative 
norm, so that each physical field has its PV counter- 
part. On the level of free Lagrangians, the physical 
and PV fields are independent from each other, while 
they are mixed by the interaction. The PV fermion 
and PV boson parts of the full Lagrangian arc 

Cb,pv = [9.$pyc>''$pv-/i?$py] , (7b) 

with mi and /ii being the PV fermion and PV boson 
masses. Note that the Lagrangians ([7]) differ by a 
minus sign from the Lagrangians (pa|) and (j2bp for the 
physical fields. The interaction Lagrangian involves 
all types of fields and has the form 

Cfb^pv ^ .go*'*'^'' + ^mvp'^-', (8) 

where 

*' = * + *py, $' = $ + $py. (9) 

The interaction is constructed in such a way that the 
physical and PV fields come into Eq. ([8]) on equal 
grounds. This feature ensures the cancellation of 
ultra-violet divergencies. The full Lagrangian com- 
bining the physical and PV contributions is thus 

Cpv = + + ^F.PV + ^B,PV + ^FB.PV- (10) 

The Lagrangian pU]) generates the interaction Hamil- 
tonian 

H'j^^ix) = -soV^VV - 5rmi;'il;' (11) 

with -0' = -0 -|- ippY and ip' = ip + ifpv- The fields ip 
and V'py cLiid ifpv) satisfy the free Dirac (Klein- 
Gordon) equations, with the corresponding masses, 
in contrast to ^ and $ which satisfy the full Heiscn- 
berg equation. The main steps leading to Eq. (jlip are 
pointed out in Appendix [XI 

The Hamiltonian pT|) has the traditional spin struc- 
ture, except for the fact that the "elementary" fields 
ip' and ip' are the sums of the physical and PV fields. 
In other words, it does not contain any contact terms 
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specific for LFD and explicitly depending on the LF 
plane orientation. This is a great merit of the PV 
regularization scheme. 

The Lagrangian ([5]), as well as the Hamilto- 
nian (fTTj) . depends on the MC 6m and on the BCC 
go- For simplicity, we consider here the case with 
only one coupling constant to be determined, but 
our scheme is completely general and can be easily 
extended to the case where many types of interaction 
occur. 

Apart from the MC and the BCC entering the origi- 
nal Lagrangian, one may also need new counterterms, 
at the level of the LF Hamiltonian, in order to re- 
store the symmetries broken by the Fock space trunca- 
tion 

[3 or by the regularization method [8| ■ We have 
already analyzed in Ref. 0] the structure of such coun- 
terterms in CLFD, using, as examples, the Yukawa 
model and QED for the case of the two-body trun- 
cation and the standard LF regularization by means 
of transversal and longitudinal cutoffs. Due to the 
explicit covariance of CLFD, the general structure of 
such counterterms can be exhibited in terms of the 
orientation, w, of the LF plane. The simplest coun- 
terterm which one may consider is given by 



lUJ-O 



(12) 



where Z^^ is a constant and l/i{LJ-d) is the opera- 
tor Eq. (jA5[) . written in covariant notations. 
This counterterm has a structure similar to that of 
the MC and appears, in all diagrams, as a factor 
Z^mi/j/{Lu-p) on each internal fermion line (here p is 
the four-momentum assigned to the line) . Other coun- 
terterms with more involved structure may appear if 
one increases the number of Fock components, giv- 
ing rise to many-body vertex corrections. The general 
renormalization scheme we propose in this paper can 
easily embrace all types of counterterms. 



the explicit form of dynamical operators, but not the 
general results discussed in this section. PV fields 
can be easily incorporated later, when we shall study 
particular physical systems. 

We are interested in the state vector, 0^'^(p), of a 
bound system. It corresponds to definite values for the 
mass M, the four- momentum p, and the total angular 
momentum J with projection a onto the z axis in the 
rest frame, i.e., the state vector forms a representation 
of the Poincarc group. This means that it satisfies the 
following eigenstate equations: 
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dp) = Pp ^'fip), (13a) 

e(p) = ^i%p), (13b) 

e(p) = -M' J{J+^) <I>'J{P): (13c) 

b^p) = M a<Pi-{p), (13d) 



where Sp is the Pauli-Lubanski vector 



c pi^ fct/3 



(14) 



and J is the four-dimensional angular momentum op- 
erator which is represented, similarly to Pp, Eq. ([5]), 
as a sum of the free and interaction parts: 

Jpii ~ Jpi) + Jpi, ■ (15) 
In terms of the interaction Hamiltonian, we have 

JpT = J H''^\x){xpuj, - XuUJp)5{Lo-x) d^x. (16) 

From the general transformation properties of both 
the state vector and the LF plane, it follows [l^ that 



JIT cl>mp) = Lp^{uj)^i^p), 



(17) 



where 



B. Covariant formulation of light-front dynamics 

In CLFD, the state vector is defined on the LF plane 
of general orientation uj-x = ^, where u; is an arbitrary 
four- vector restricted by the condition uj'^ ~ 0, and ^ is 
the LF "time". We shall take .J = 0, for convenience. 

Let us recall here, for completeness, how the state 
vector of a physical system is constructed. In order 
to avoid congesting notations, we do not consider for 
the moment PV fields. These fields infiuence only 



(18) 



The equation (jl7p is called the angular condition. We 
can now use it in order to replace the operator J™* 
entering into Eq. ([T3]) by Lp^(w). Introducing the no- 
tations 



Mp. = Jfj+Lp,{uj), 



Wp = -ep^^p Af"^ 



(19a) 
(19b) 
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we obtain, instead of Eqs. (|13cp and (jl3dp : 

W^^^ip) = -M^J{J+l)c^'j{p). (20a) 
VK3 (l^i^p) = M a cpi%p). (20b) 

These equations do not contain the interaction Hamil- 
tonian, once (f> satisfies Eqs. (|13ap and (jlSbp . The 
construction of the wave functions of states with def- 
inite total angular momentum becomes therefore a 
purely kinematical problem. Indeed, the transforma- 
tion properties of the state vector under rotations of 
the coordinate system is fully determined by its total 
angular momentum, while the dynamical part of the 
latter is separated out by means of the angular con- 
dition. The dynamical dependence of the wave func- 
tions on the LF plane orientation now turns into their 
explicit dependence on the four- vector to Q. Such 
a separation, in a covariant way, of kinematical and 
dynamical transformations is a definite advantage of 



where 4'n,ija'{- ■ ■) is the n-body LF wave function 
(Fock component) describing the state made of one 
free fermion and (n— 1) free bosons, a'^ {c^) are the free 
fermion (boson) creation operators, e^, = -^/k^ + m^, 
and mi is the mass of the particle I with the four- 
momentum ki. The combinatorial factor l/(n— 1)! is 
introduced in order to take into account the identity 
of bosons.^ The variables t„ describe how far off the 
energy shell the constituents are. As explained in Ap- 
pendixlBl the momentum lotu can be identified with a 
fictitious particle, called spurion. In practical calcula- 
tions, the infinite sum over n is truncated by retaining 
terms with n which does not exceed a given number 



^ Usually the factor I/a^/ (n — 1)! is used, instead of l/(n — 1)!. 
Our choice however allows to remove additional combinato- 
rial factors in the equations for the Fock components, which 
would arise in the former case. 



CLFD as compared to standard LFD on the plane 
t + z = 0. 



C. General Fock decomposition of the state 
vector 

According to the general properties of LFD, men- 
tioned in the Introduction, we decompose the state 
vector of a physical system in Fock sectors. Schemat- 
ically, we have 

^i%p)^\l) + \2) + --- + \n)^... (21) 

Each term on the r.-h.s. denotes a state with a fixed 
number of particles from which the physical system 
can be constructed. In the Yukawa model the analyt- 
ical form of the Fock decomposition is 



(22) 



N , while those with n > N are neglected. Decompo- 
sitions analogous to Eq. can be easily written for 
the QED case 0| or for a purely scalar system [g^l . 

The normalization condition for the state vector is 
given by 



= 2ep 6,^,,S^^Hp - P')- (23) 



Being rewritten through the CLFD wave functions, it 
has the form 



00 



Op) = / 0„...'(fci...A:„,p,c.r„)at,(ki)ct(k2)...ct(k„)|O) 



S^"^^ (ki + . . . + kn - p - UJTn)2{uj-p)dTn 



d^k. 



\ (2^)3/2^2^' 
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where 

In = 



(uj-p) 



(27r)3("-i)(n- 1)! 



n.(T(j Tu.aa 



(25) 



is the relative contribution of the 7i-body sector to the 
full norm. For shortness, we omitted the arguments of 
the wave functions. The factor 1)! in Eq. 

appears as a combined effect caused by the presence 
of the same factor in Eq. and by the contrac- 
tion of the creation and annihilation operators, when 
calculating the l.-h.s. of Eq. 

As follows from the discussion in Sec. Ill Bl the spin 
structure of the wave functions 4>n,aa' is very simple, 
since it is purely kinematical, but it should incorpo- 
rate w-dcpendcnt components in order to fulfill the 
angular condition p7p . It is convenient to decompose 
each wave function ^n.o-o-' into invariant amplitudes 
constructed from the particle four-momenta (includ- 
ing the four-vector uj\) and spin structures (matrices, 
bispinors, etc.). In the Yukawa model, for instance, 
we have 



t>2,<y<y' = Wdr'(fcl) 



-02 +-02 



U}-p 



(26a) 
Ua[p), (26b) 



since no other independent spin structures can be con- 
structed. Here it's are bispinors, ipi, -02, and ?/'2 are 
scalar functions determined by the dynamics. For a 
spin 1/2 system coupled to scalar particles, the num- 
ber of invariant amplitudes for the two-body Fock 
component coincides with the number of indepen- 
dent amplitudes of the reaction spin 1/2-1- scalar 
spin 1/2-1- scalar, which is (2 x 2)/2 = 2, due to parity 
conservation. 



Note that the formulas (|22l), and ^ are writ- 
ten for the state vector which contains physical par- 
ticles only. The use of the PV regularization, strictly 
speaking, changes them. However, their generaliza- 
tion is straightforward. We do not give here the cor- 
responding general relations, but give their particular 
forms when we proceed to the consideration of con- 
crete physical systems. 



D. Eigenstate equation 

The equations for the Fock components can be ob- 
tained from Eq. (jl3bp by substituting there the Fock 
decomposition (|22[) of the state vector 0(p) (here and 
below wc will omit, for shortness, all indices in the no- 
tation of the state vector) and calculating the matrix 
elements of the operator in Fock space. With the 
expressions ([5]) and (O, we can easily get the eigen- 
state equation [20| : 



2{u:-p) / ff™*(c.r)^0(p) = - 



(27) 

where i^*"* is the interaction Hamiltonian in momen- 
tum space: 



(28) 



For the Yukawa model with the PV regularization. 
i7™*(a:) is given by Eq. (fTTI) . 

According to the decomposition ([2^ . the conserva- 
tion law for the momenta in each Fock component has 
the form 



ki+ h A;„ = p + WT„ 



(29) 



Hence, the action of the operator — AI on the 

state vector reduces to the multiplication of each Fock 
component by the factor — = 2(u;-p)r„. 

It is therefore convenient to introduce the notation 



g{p) = 2{uj-p)t(I){p), 



(30) 



where f is the operator which, acting on a given com- 
ponent 4)n,aa' of 0(p), givcs r„0„^crcr' . G (p) has the 
Fock decomposition which is obtained from Eq. ([22|) 
by the replacement of the wave functions (pn^aa' by 
the vertex functions r„ (which we will also refer to as 
the Fock components) defined by 



(31) 



and s„ = (fci + . . . knY ■ Since for each Fock compo- 
nent s„ — = 2{iu-p)Tn, we can cast the eigenstate 
equation in the form 



Zn 



-W^\lot) —g{p) . (32) 



The physical bound state mass M is found from the 
condition that the eigenvalue is 1. This equation is 
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quite general and equivalent to the eigenstate equa- 
tion (fT3b|) . It is non-perturbative. 

The normalization integrals (j25p rewritten through 
the vertex functions are 

(27r)3("-i)(n- 1)! 

with f „ = 7°rt 7°. 

Since the Hamiltonian ()lip has the ordinary struc- 
ture in terms of the fields ip' and ip' (i.e. it does not 
include any contact terms), we may apply, for sub- 
sequent calculations, the standard CLFD graph tech- 
nique rules @ with minor changes (see Appendix iBjl . 



III. SYSTEMATIC RENORMALIZATION 
SCHEME IN CLFD 

In the usual renormalization scheme, the bare pa- 
rameters^ are determined by fixing some physical 
quantities like the particle masses and the physical 
coupling constant. The physical parameters are thus 
expressed through the bare ones. This identification 
implies in fact the following two important conse- 
quences which are usually never clarified in LFD cal- 
culations, but are at the heart of our scheme. 

(i) In order to express the physical parameters 
through the bare ones, and vice versa, one should be 
able to calculate observables or, in other words, phys- 
ical amplitudes. In LFD, any physical amplitude is 
represented as a sum of partial contributions, each 
depending on the LF plane orientation. Since an ob- 
servable quantity can not depend on the choice of the 
LF plane, this spurious dependence must cancel in the 
whole sum. Such a situation indeed takes place, for 
instance, in perturbation theory, provided the regular- 
ization of divergencies in LFD amplitudes is done in 
a rotationally invariant way Q. In non-perturbative 
LFD calculations which are always approximate (say, 
due to the Fock space truncation we just use here) 



^ The term " bare parameters" means here the whole set of pa- 
rameters entering into the interaction Hamiltonian, e.g. the 
BCC, the fermion MC, etc. 



the dependence on the LF plane orientation may sur- 
vive even in physical amplitudes. For this reason, 
the identification of such amplitudes with observable 
quantities becomes ambiguous and expressing the am- 
plitudes through the physical parameters turns into a 
non-trivial problem. 

When working in standard LFD on the plane t+z = 
0, one may think that even approximate LFD ampli- 
tudes do not depend on the LF plane orientation. As a 
matter of fact, they do, but this dependence is simply 
hidden. It reflects itself by the non-invariance of the 
corresponding amplitudes: the result of calculation is 
affected by the choice of the reference frame. 

As we have detailed in Sec. [TTl CLFD is a unique 
tool to control this dependence in terms of the arbi- 
trary light-like four-vector uj specifying the LF plane 
orientation. We shall see in Sec. IIVI how one should 
make use of this property to define the physical 
fermion-boson coupling constant from the two-body 
Fock component of the type (|26b[) . 

(ii) The explicit form of the relationship between 
the bare and physical parameters depends on the ap- 
proximation which is made. This is trivial in pertur- 
bation theory where the order of approximation is dis- 
tinctly determined by the power of the coupling con- 
stant. In our non-perturbative approach based on the 
truncated Fock decomposition an analogous parame- 
ter is absent. At the same time, to make calculations 
compatible with the order of truncation, one has to 
trace somehow the level of approximation. This im- 
plies that, on general grounds, the bare parameters 
should depend on the Fock sector in which they are 
considered. Moreover, this dependence must be such 
that all divergent contributions are cancelled. 

We will show in the following how to realize the 
renormalization procedure in practice. For clarity, we 
take, as a background, a model of interacting fermions 
and bosons like the Yukawa model or QED. At the 
same time, let us emphasize once more that the gen- 
eral renormalization strategy developed here is appli- 
cable to physical systems with arbitrary interaction 
admitting Fock decomposition of the state vector. 

A. Fock sector dependent counterterms 

Let us call the maximal number of Fock sectors 
considered in a given approximation, and n the num- 
ber of constituents in a given Fock sector [one fermion 
and (n—l) bosons]. We have n < N. Each Fock sector 
is described by its vertex ftmction defined by Eq. ([?T|) . 
In a truncated Fock decomposition, each vertex func- 
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tion should depend on N. We shall thus denote it 
by r^^'' . Graphically, rl^'' can be represented by the 
diagram shown in Fig. [31 




(n — 1) bosons 



FIG. 3: Vertex function of order n for the A'^-body Fock 
space truncation. 



1. Mass counterterm 

The simple example of the fcrmion self-energy 
renormalization by the MC within the two-body Fock 
space truncation, presented in the Introduction, can 
serve as a guideline to define our general rules. In 
this example, the MC should be labelled with a sub- 
script and denoted by dm2, in order to indicate that 
it is introduced in order to cancel, at ^ = to, the 
fermion self-energy contribution which belongs to the 
two-body Fock sector. In other words, Sm2 is related, 
by construction, to the two-body state, even though 
it is attached to a single fermion line. More gener- 
ally, the subscript at 5m corresponds to the maximal 
number of particles in which the fermion line where 
the MC is attached can fluctuate, so that the total 
number of particles at any LF time equals A^. In the 
given example it is 2. 

Let thus denote by 5mi the MC in the most general 
case. Since we truncate our Fock space to order TV, 
one should make sure that, at any LF time, the total 
number of particles is at most N. Our first rule is 
thus: 

• in any amplitude where the MC 5mi appears, 
the value of I is such that the total number of 
bosons in flight plus I equals the maximal number 
of Fock sectors considered in the calculation, i.e. 
N. 

For instance, in the typical contribution indicated 
in Fig. 131 the MC is (5m(^_„+i). Indeed, since there 
are already (n — 1) bosons in flight, the fermion line 
can fluctuate in at most I = {N — n + 1) particles, so 



that the total number of particles at a given LF time 
is just {N — n + 1) + {n — 1) = N and no more. 




(n — 1) bosons 



FIG. 4: Typical MC insertion. 

In order to calculate the whole set of the MC's Smi 
one should proceed in the following way. Any cal- 
culation of the state vector within the Fock space 
truncation of order N involves the MC's Smi with 
I = 1,2, ... N. We emphasize at this point that the 
MC's 6mi are successive approximations to the true 
MC Sm appearing in the original Hamiltonian. They 
are connected to each other by some kind of recur- 
sion in the sense that finding 6mi requires knowledge 
of the lower order counterterms, i.e. Smi, 5m2, .. 
Smi^i. This is analoguous to what happens in any 
perturbative calculation where each MC relates to a 
definite order of perturbation theory. 

For the MC of lowest order, we simply have 

5mi = 0, (34) 

which is trivial since if the fermion can not fluctuate 
in more than one particle, its mass is not renormal- 
ized at all. The subsequent MC's are calculated by 
solving successively the eigenvalue equations for the 
Fock components for the 2,3, . . . A^-body Fock space 
truncations. 



2. Bare coupling constant 

Let us now come to the determination of the BCC. 
The general strategy we developed above for the cal- 
culation of the MC should of course be also applied to 
the BCC case, with however a bit of caution, since the 
BCC may enter in two different types of contributions. 

• The first one appears in the calculation of the 
state vector itself, when Eq. ([32]) is solved. In 
that case, any boson-fermion coupling constant 
is associated with the emission or the absorp- 
tion of a boson which participates in the particle 
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counting, in accordance with the rules detailed 
above, since it is a part of the state vector. 

• The second one appears in the calculation of 
the boson-fermion scattering amplitude or of 
the boson-fermion three-point Green's function 
(3PGF) like the electromagnetic form factor in 
the case of QED. This observable is usually con- 
sidered, at some kinematical point, to define the 
physical coupling constant. Now the external 
boson is an (asymptotic) free field rather than 
a part of the state vector. The particle count- 
ing rule advocated above should therefore not 
include the external boson line. 

One has thus to distinguish two types of BCC's: 
go and go describing, respectively, the interaction ver- 
tices of the fermion with internal and external bosons. 
As we shall see below, these two BCC's are found 
from different conditions and do not coincide with 
each other for a finite order Fock space truncation. 
The necessity to introduce two BCC's can be also cx- 



e e 




/ / 

FIG. 5: Typical contribution to the two-body scattering 
amplitude of a probe on a bound state system. 

plained from another point of view. Let us consider 
the scattering amplitude of a given probe on a bound 
state system. Such scattering amplitude can be rep- 
resented by the diagram in Fig. The in and out 
asymptotic states are defined in terms of the (struc- 
tureless) probe denoted by e and the bound state sys- 
tem denoted by /. The state vector of the bound sys- 
tem is calculated within a given approximation (the 
Fock space truncation, in our case), starting from a 
known Hamiltonian. Therefore, the calculated state 
vector "knows" nothing about the subsequent inter- 
action it can have with the probe, i.e. it should be 
independent of any coupling to the external virtual 
bosons exchanged between the probe and the bound 
system. 

Similarly to the MC, the BCC's should also keep 
track of the Fock sector in which they appear. To il- 
lustrate this fact, let us write down some typical con- 
tributions to the fermion self-energy, involving at most 




FIG. 6: Two- and three-body Fock sector contributions to 
the fermion self-energy. 



two bosons in flight (i.e. for iV = 3). They are shown 
in Fig. [51 All the vertices are described by the in- 
ternal BCC's, since the self-energy is a part of the 
fermion state vector. The vertices (1) and (2) involve 
the BCC's attached to the two-body sector in three- 
body truncated Fock space. So, each of these vertices 
or both can be "dressed" by one more bosonic loop, 
as indicated for the vertices (4) and (5). The latter 
vertices correspond to states fully "saturated" with 
bosons. In other words, no radiative corrections to 
them arc allowed in the given approximation. From 
here it follows that the vertices (l)-(3) and (6), on 
the one hand, and (4) and (5), on the other hand, 
are described by different BCC's. Analogously to the 
MC, we will denote each internal BCC by goi, where 
the subscript I indicates which Fock sector the given 
BCC belongs to. We can then formulate the general 
rule: 

• in any amplitude which couples constituents in- 
side the state vector one should attach to each 
vertex the internal BCC goi. The value of I is 
such that the total number of bosons in flight be- 
fore (after) the vertex - if the latter corresponds 
to the boson emission (absorbtion) - plus I equals 
the maximal number of the Fock sectors consid- 
ered in the calculation, i.e. N. 

Applying this rule to the diagrams shown in Fig. [6l 
we find that the vertices (l)-(3) and (6) are described 
by the BCC goa [no bosons in flight before the ver- 
tices (1) and (3) or after the vertices (2) and (6)], 
while the vertices (4) and (5) are described by the 
BCC ,902 (one boson in flight before or after each ver- 
tex). The BCC go2 is calculated within the two-body 
Fock space truncation and enters into the calculations 
within the three-body truncation as a known parame- 
ter, while go3 should be found. When considering the 
four-body truncation, one should then calculate go4, 
knowing both BCC's 502 and gos. The procedure can 
thus be extended to arbitrary A^. 

Let us finally consider typical contributions to the 
3PGF, analogous to those for the self-energy, as dis- 
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FIG. 7: Typical contribution to the fermion-boson 3PGF. 



cussed above. Some of them are presented in Fig. [T] 
The crucial difference with the self-energy case con- 
sists in that the bosons emitted from the vertices (1)- 
(3) are absorbed by an external particle (typically an 
external probe) which is not included into the state 
vector. For this reason, the vertices (l)-(3) are de- 
scribed by the external BCC's goi, while all the others 
correspond to the internal BCC's goi. The values of 
/ for the "internal" vertices are determined according 
to the rule specified above, with the only exception 
that the external boson does not take part in the par- 
ticle counting. For the "external" vertices (l)-(3) the 
situation is quite similar. The vertices (1) and (2) cor- 
responding to the two-body Fock sector (the external 
boson is not counted!) admit additional corrections, 
within the three-body Fock space truncation, from 
"internal" bosonic loops, while the vertex (3) does not, 
since it appears within the three-body sector. We can 
conclude from here that the external BCC's attached 
to the vertices (1) and (2) should differ from that for 
the vertex (3). We can thus formulate the following 
general rule: 

• in any amplitude which couples constituents of 
the state vector with an external field, one should 
attach to the vertex involving this external field 
the external BCC goi- The value of I is such 
that, at the LF time corresponding to the ver- 
tex, the total number of internal bosons in flight 
- i.e. those emitted and absorbed by particles 
entering the state vector - plus I equals the max- 
imal number of the Fock sectors considered in 
the calculation, i.e. N. 

This rule prescribes to attach to the vertices (1) and 
(2) in Fig. [7] the BCC 502, while the vertex (3) is 
associated with 501 • 

The lowest order BCC's are 

501 = 0, (35a) 
501 = 5- (35b) 

Eq. (|35ap is trivial, because no fermion-boson interac- 
tion is allowed in the one-body Fock space truncation. 



Eq. (j35bp reflects the fact that the external BCC, in 
the same approximation, is not renormalizcd at all 
since a single fermion can not be " dressed" . 

Note that the rule for attaching BCC's to exter- 
nal vertices holds unchanged if the external field is 
a particle of another sort than the bosons entering 
into the state vector. Such a situation takes place, for 
instance, when one calculates the electromagnetic ver- 
tex of a fermion, when the state vector does not con- 
tain photons (e.g. in the Yukawa model). Evidently 
no problems with particle counting arise in this case. 

Some illustrations of the rules concerning the inter- 
nal and external BCC's are given in Fig. [S] 




(n — 1) bosons 



FIG. 8: Typical contributions to the fermion state vector 
for the absorption (a) and the emission (b) of an internal 
boson, and to the fermion-boson 3PGF (c). 



This completes the set of our general rules to define 
in a systematic and non-perturbative way the MC and 
the BCC's in LFD calculations within truncated Fock 
space. The three rules exposed above have very sim- 
ilar logical grounds based on the particle counting in 
intermediate states. Namely, the index I at the MC 
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and the BCC's is always calculated by the same rule: 
I equals the difference between the order of approxi- 
mation N and the number of internal bosons in flight 
at the corresponding LF time. 

Though we rely on the fermion-boson model when 
weeonsidered the above procedure, this latter can be 
easily extended to other systems with additional coun- 
terterms and bare parameters. 

B. Renormalization Conditions 

Once proper bare parameters and counterterms 
have been identified, one should fix them from a set 
of renormalization conditions. In perturbation the- 
ory, there are three quantities to be determined: the 
fermionic MC, the BCC, and the norm of the fermion 
field (if the boson field renormalization is not con- 
sidered). Usually, the on-mass-shell renormalization 
is applied, with the following conditions. The MC is 
fixed from the requirement that the two-point Green's 
function has a pole at = M^, where M is the phys- 
ical mass of the fermion. The fermion field normal- 
ization is fixed by the condition that the residue of 
the two-point Green's function in the pole is 1. The 
BCC is determined by requiring that the on-mass-shcll 
3PGF is given by the product of the physical coupling 
constant and the elementary (i.e. not "dressed") ver- 
tex. 

One should now extend these conditions in order to 
determine the bare parameters and the counterterms 
in a non-perturbative LED framework. We do not 
need to rcnormalize the fields in our approach, since 
we deal with already normalized state vectors. It re- 
mains to specify the procedure to find the set of MC's 
Smi and the BCC's gai and goi. 

The eigenstate equation for the state vector in the 
TV-body approximation includes two unknown param- 
eters, the MC Stun and the BCC goN, while all Smi 
and go; with Z < iV — 1 are defined from calculations 
made for Fock space truncations of lower orders. Pro- 
vided goN is fixed, the MC Sttin is foimd from the con- 
dition that the eigenvalue in Eq. ([5^ equals 1 in the 
limit where the mass of the ground state, M, is equal 
to the mass of the fermionic constituent, m. To de- 
termine goN, one should relate it to the physical cou- 
pling constant. Eor instance, in the traditional QED 
renormalization scheme, the BCC is found from the 
requirement that the residue of the photon-electron 
scattering amplitude in the pole s = m? , where s is 
the invariant energy squared of the system electron + 
photon, would be proportional to the physical elec- 



tron charge squared. This condition is not restricted 
to perturbation theory and can be directly extended 
to non-perturbative approaches as well. In our lan- 
guage, it is equivalent to require that the two-body 
vertex F2 is proportional, at s = m^, to the elemen- 
tary vertex, the proportionality coefficient being just 
the physical coupling constant. Eor example, one has 
to require, at s = m^, that E2 = g in the Yukawa 
model, and F2 = 67" in QED. Erom here a relation 
between the internal BCC and the physical coupling 
constant can be derived. 

To determine the external BCC, a similar condition 
should be imposed on the 3PGE. Indeed, the two-body 
vertex and the 3PGE represent two different channels 
of the same reaction. Hence, taken entirely on the 
mass and energy shells, they must coincide. It follows 
from here that the 3PGE G{p,p'), where p [p') is the 
four-momentum of the incoming (outgoing) fermion, 
reduces, at p^ = p' = and = {p' — p)^ = n^, 
to the same product of the elementary vertex and the 
physical coupling constant. Since the internal BCC 
has already been found from the calculation of 
the above condition on the 3PGE allows to relate the 
external BCC with the physical coupling constant. 

Note that the analytical continuation of the two- 
body vertex function and the 3PGE to the non- 
physical points s = m? and = /i^ does not en- 
counter any technical difficulties, even in numerical 
calculations, since we can use for this aim the eigen- 
state equation ([5^ . Its l.-h.s. contains the function 
to be continued in the non-physical point. Erom the 
r.-h.s., both F2 and 3PGF are expressed through in- 
tegrals involving the vertex functions in the physical 
domain, whereas the dependence of the integrand on 
s and q^ is explicit. We can put there s ~ m? and 

In perturbation theory, the equivalence of the on- 
mass-shell two-body vertex function and 3PGE (cal- 
culated in the same order!) appears automatically, as 
a consequence of the analytical properties of scatter- 
ing amplitudes. Eor this reason, it does not make any 
sense to distinguish the external and internal BCC's, 
because they are equal to each other to any order. The 
same would happen in exact calculations, if they were 
possible. In our non-perturbative approach based on 
truncated Eock decompositions, the non-renormalized 
two-body vertex function and the 3PGE, even taken 
on the energy shell, do not automatically coincide , in 
any finite order approximation. Moreover, they may 
be not constant (i.e. keep dependence on particle mo- 
menta) and depend also on the LE plane orientation 
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(on the four- vector w, in CLFD). If so, the question 
how to identify such objects with physical constants 
in the renormaUzation point requires special consider- 
ation. 

First of all, one has to fix unambiguously all kine- 
matical variables in the renormalization point, in or- 
der that both the two-body vertex function and the 
3PGF would turn into constants. Concerning the 
3PGF, we can offer a universal solution of this prob- 
lem. Indeed, G{p,p') depends dynamically on two 
scalar variables, say, and to-q/uj-p. The condition 
q"^ = ^j? fixes the first variable, but leaves free the 
second one. However, if one imposes the condition 
Lo-q = (analogous to = in ordinary LFD) on 
the four-vector w, the dependence of the 3PGF on the 
second variable drops out, and it becomes a constant 
at fixed q^. This condition is also necessary for the fac- 
torization of the total scattering amplitude of a probe 
on the physical system under study, in terms of the 
external boson propagator and the 3PGF. In practice, 
one should calculate G{p,p'), keeping ui-p = uj-p' , and 
then continue it analytically, as a function of only one 
variable to the point q^ = A*^- 

The two-body vertex r2 also depends dynamically 
on two scalar variables, e.g. s and x = uj-k2/uj-p, 
where ^2 is the boson four-momentum. The condition 
s = TV? does not fix x. In an exact calculation, r2 at 
s = m? does not depend on x, while the x-dependence 
survives because of approximations. Therefore, some 
additional restriction must be imposed on x. Unfortu- 
nately, a universal choice how to fix x is hardly possi- 
ble, in contrast to the case of the 3PGF. This problem 
should be solved separately, for each particular phys- 
ical system. 

Once the kinematics in the renormalization point 
is fixed, both the two-body vertex function and the 
3PGF turn, in this point, into constant matrices. 
Their dependence on the LF plane orientation may 
however survive. For example, in the Yukawa model, 
the term "02 "t- (^/(^'P) in the two-body wave func- 
tion (j26bp . which explicitly depends on implies 
analogous dependence in the two-body vertex func- 
tion, unless ^2 = ill ftic renormalization point. To 
get rid of it, one should insert new countcrtcrms into 
the LF Hamiltonian, also explicitly depending on w, 
which cancel completely the cj-dependent term in r2. 
Additional counterterms may also be needed to kill 
possible w-dependence of the 3PGF. According to the 
general renormalization strategy, these counterterms 
must depend on the Fock sector under consideration, 
in full analogy with the MC and the BCC's. 

Introducing the internal and external BCC's and 



imposing the above renormalization conditions on 
both the on-energy-shell two-body vertex and 3PGF, 
we force their coincidence, after renormalization, for 
arbitrary Fock space truncation of finite order. At 
this level, we restore the cross-invariance of scattering 
amplitudes in the renormalization point. 

Summarizing, we propose the following non- 
perturbative renormalization conditions: 

• The MC is fixed by solving the eigenstate equa- 
tion h32j) in the limit AI m, where M is the 
mass of the ground state of the physical system, 
and m is the physical mass of the fermionic con- 
stituent. 

• The state vector is normalized according to the 
standard condition i23\) . 

• The internal bare coupling constant g^i is fixed 
from the condition that the uj -independent part 
of the two-body vertex function taken at s = n? 
and at a given value of x, denoted by x* , is 
proportional to the elementary vertex, with the 
proportionality coefficient being the physical cou- 
pling constant. 

• The external bare coupling constant goi is fixed 
from the condition that the uj -independent part 
of the 3PGF calculated at cu-q = and taken at 
q^ = fj? is proportional to the elementary ver- 
tex, with the proportionality coefficient being the 
physical coupling constant. 

• The uj-dependent counterterms in the Hamil- 
tonian are fixed by the conditions that the u- 
dependent parts of the two-body vertex function 
in the point (s = m'^,x*) and the 3PGF in the 
point {q^ = pi?,ijj-q = 0) are zero. 

• The values of all bare parameters and countert- 
erms for I < N are determined from succes- 
sive calculations within the 1, 2, ... TV Fock space 
truncations. 



IV. APPLICATIONS TO THE N=2 FOCK 
SPACE TRUNCATION 

In order to show simple but nevertheless meaningful 
applications of the general renormalization scheme de- 
veloped above, we consider, in a first step, the Yukawa 
model. The LF Hamiltonian including PV fields has 
been derived in Sec. [Hi We shall then extend our 
results to QED which is very similar to the Yukawa 
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model for the N = 2 Fock space truncation. To sim- 
plify notations, we will omit in the following the su- 
perscript (N) at each vertex function. 



Yukawa model 



Solution of the eigenstate equation 



The solution of the eigenstate equation p2l) within 
the N ~ 2 Fock space truncation has already been 
found in Rcf. Q- Sharp cut-offs imposed on the lon- 
gitudinal and transverse components of particle mo- 
menta were used to regularize the amplitudes. We 
revisit here the same problem, but apply the PV reg- 
ularization method (with one PV fermion and one PV 
boson), as advocated in SecHIl Besides that, we shall 
apply and test the new rcnormalization scheme pro- 
posed in Sec. IIIII 

The use of the PV rcgularization extends Fock 
space: instead of one Fock component for the one- 
body sector we have now two components correspond- 
ing to the physical and PV fcrmions, while the two- 
body sector is described by four components related 
to the states either with both physical particles, or 
with the physical fermion plus the PV boson, or with 
the PV fermion plus the physical boson, or, finally, 
with both PV particles. The extension of Fock space 
makes the computational analysis more cumbersome, 
but this is compensated by the simplification of the 
equations for the Fock components due to the disap- 
pearance of the contact terms, as well as the absence 
of spurious w-dcpendent contributions to the fermion 
self-energy Q. 

To incorporate the PV particles into the state vec- 
tor, we will supply the vertex functions, as well as the 
particle momenta and masses, with the indices i and 
j which relate, respectively, to fermions and bosons. 
Each index is either for a physical particle or 1 for 
a PV one. All particle momenta are on their mass 
shells: 

1,2 _ 2 ; 2 _ 2 

with niQ — m and fiQ = fi being the physical particle 
masses. 

Following Eq. (j26bp . we decompose the vertex func- 
tions in invariant amplitudes: 

uipu)n<p) = {m\ - M^)^Plu{pu)uip), (36a) 



u{ki 



u{ku) 



mi/) 
uj-p 



u{p). (36b) 



Since the fermion momenta for the one- and two-body 
vertices are different, we denote them by different let- 
ters, namely, by pi, for the one-body vertex and by kn 
for the two-body one. The boson four-momentum is 
k2j ■ Here M is a temporary notation for the physical 
fermion mass (in the end of the calculation we will take 



the limit M — > m) 



k"^ 



ipl are constants 



(i.e. they do not depend on particle momenta), while 
the invariant functions and 62"' may have momen- 
tum dependence.'^ We introduce the standard LFD 
variables 



2j 



LU-p 



R± = k2j± - xp 



(37) 



which are, respectively, the longitudinal and trans- 
verse (with respect to the LF plane orientation) com- 
ponents of the bosonic momentum. Note that the 
square of the two-dimensional vector is an invari- 
ant: i?^ = — (fc2j — xp)'^. 

The system of coupled equations for the vertex func- 
tions in the two-body truncated Fock space is shown 
graphically in Fig. [HI Since the Hamiltonian pT|) does 



5m2 



902 




S02 



FIG. 9: System of equations for the vertex functions in 
the two-body approximation. 



not include contact terms, this system is much sim- 
pler than its analog considered in Ref. For sim- 
plicity, we have not drawn the lines associated with 
the spurion (see discussion in Appendix IB|) . Besides 
that, we do not need to introduce any specific coun- 
tertcrms which explicitly depend on the LF plane ori- 
entation. These are evident merits of the PV rcgu- 
larization. The amplitudes of the LF diagrams are 
calculated according to the graph technique rules (see 
Ref. @ and Appendix [B|) . We denote the intermedi- 
ate fermion four-momenta in the one- and two-body 



^ In the two-body approximation, as wc will see, they reduce 
to constants, but this is no more the case in higher order 
approximations. 
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states as p[^, and k[^,, respectively. The intermedi- where 
ate boson four-momentum is fcjj' ■ The values of the 
intermediate particle momenta are defined from the 
conservation laws in the vertices, taken for n = 1 
and n = 2. 

The system of equations reads 



u{pu)r\u{p) = uipu){Vi+V2)u{p), (38a) 
uiku)r^'uip) = uiku)V3uip), (38b) 



where 



dx' 



(27r)3 7o 2x'{l-x') 



E 



(-1) 



'12 



(41) 



which is nothing else than the PV-regularizcd fermion 
two-body self-energy (apart from the coupling con- 
stant). After integration, it depends on the four- 
momentum p only. We will use the following decom- 
position: 



Vi 



bm-. 



502 



M2 

rf^i?' /•! dec' 



(27r)3 7o 2x-'(l-x') 

502 2^(-l) _2 _ ,.2^1, 



(39a) 
(39b) 



(39c) 



and s^l 



(fcj^j, + fc2j')^- We took into account that 
(5mi = 0. Necessary kinematical relations used for 
deriving Eqs. p9p are given in Appendix I C II We in- 
troduced the integration variables R'^ and x' related 
to the intermediate boson momentum k'^y in full anal- 
ogy with Eqs. (|37|) . The summations in Eqs. ((39)) run 
through to 1 in each index. The system of equa- 
tions (l38l) must be solved in the limit M m. 



We have four unknown functions b^^ , four unknown 
functions b2 and two unknown constants ip{. So, we 
have to deal with a system of 10 x 10 linear integral 
equations. However, Eqs. (|38bp and (|39c|) allow to ex- 
press easily Fj'' through T\. Since T\ are constants, 
it follows from these equations that bi and &2"' 
constants too and they depend neither on i, nor on j. 
Then, due to the fact that the vertex functions come 
into all equations, being sandwiched with on-mass- 
shell bispinors, we may simply substitute every- 



where by the quantity V3, Eq. (|39cp . It thus follows 
that we can rewrite Eq. (|39b[) in the form, in the limit 
M ^ m 



V2 



'902 S(p)| 



=M'- 



v. 



(40) 



S(p) s Aip') 



m 



(42) 



The explicit form of the functions A and B for arbi- 
trary values of p^ can be found in Ref . . 

Now, substituting Eq. (|39cp into Eq. (gD]) and then 
into Eq. p8a|) . and using that r\ = {mf - is 
proportional to the unity matrix, we turn the latter 
equation into a system of two linear matrix equations 
for ipi and ipl . Multiplying these equations by u{p) to 
the right and by u{pii) to the left, summarizing over 
the fermion polarizations, and taking the trace, we ar- 
rive at a system of two linear homogeneous equations 
for the one-body vertices. This system reads 



= 



(43) 



ith 



^00 



,11 



-2(m + mif {5m2 - .go2(^ + B)} , 
2(m + mif {(57712 - 5o2(^ + B)} , 

- |m — mmi (mi -f 2om2) 



m 
2 



\ 2Animi + B{m 



where A — A{m^), B — B{m^). Equating the deter- 
minant of the system to zero, we obtain a quadratic 
equation for 5m2 with the solution 



5m2^9l2{A + B). 



(44) 



The second root is rejected because it does not disap- 
pear at (702 = 0. Substituting Eq. ([44]) into any of the 
two equations (|43|) we find 



ai. 



= 0, 



(45) 
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where ai is an arbitrary constant. Substituting this 
solution into Eq. (|39cp . then into Eq. (j38bp . and com- 
paring the resuh with the r.-h.s. of Eq. (j36bp . we 
easily obtain 



0. 



(46) 



It is interesting to compare the solution (|l5|) 
and for the one- and two-body vertex functions 
with that found in Ref. Q- Due to the extension of 
the Fock space basis, the vertex functions considered 
here have more components, but both solutions pos- 
sess nevertheless the same main features. Firstly, they 
are constants, i.e. momentum independent. Secondly, 
the one-body vertex has only one component: the one- 
body wave function of the PV fermion, -0}^, vanishes 
identically, while the physical component is not 
fixed and must be computed from the normalization 
condition for the state vector. Thirdly, the two-body 
vertex has no w-dependent part, since the components 
are zero. Finally, the form of the solution (|46p is 
the same as that from Ref. 0| , apart from the coupling 
constant which is now go2 instead of g. The same is 
true for the MC 



2. Normalization of the state vector 

The Fock components found above must be prop- 
erly normalized. As we have already explained, the 
formulas and ([55)) for the " partial" normalization 
integrals must be modified, since we have to take into 
account the sectors which contain PV particles. This 
is done very easily. One should simply sum over all 
possible two-body states, keeping in mind that each 
PV particle brings the factor (—1) to the norm. The 
contribution of the one-body state is thus 

^ ^( ;^)» Tr{(^ + m)(^i, + m,)} ^^^^s^2 



4m2 {i,iy^Am\l 



(47) 



We used here Eqs. (jClap and (|C2p at A/ — m from 
Appendix IC II The norm of the two-body state reads 



(6i^yTr{(^ + TO)(^i,+m,)} 



i+j 



(sj2 — m?)^ 



(48) 



where 



1 -a; 



Substituting 6f = 277150211 into Eq. ([^5]) . calculating 
the trace and using Eqs. (|Clc[) . (jC3p . and (|C5|) . we 
obtain 



R'i + [(1 - x)m + m^f 



[R\ + (1 — x)^'^A + xmf — x{l — x)m? 



(49) 



Note that in order to get a finite result for I2, the 
bosonic PV regularization is enough. For this reason, 
one can retain in the sum over i the term with i = 
only^ : 



I2 = 4m^ai5o2 J2, 



(50) 



where 



J2 = 



87r2 



dR I R I 



dx i 



E 



{-iy[Rl + {2-xY 



[R\ + (1 - .T)Ai| + x^m^ 



(51) 



From the normalization condition (j24[) in the two- 
body approximation one gets Ii + 12 = \ and, hence, 

^1 = A 2n ! 2 M - (52) 

Am^[l+g^^J2) 

The integral J2 diverges logarithmically when the PV 
boson mass /ii tends to infinity: 



J2 = ^l0g^ + /(ii), 
lOTT"' TO Vm/ 



(53) 



where / is a finite function independent of ^i. Terms 
vanishing in the limit /ii — > 00 are neglected. Note 
that if /.Ji is large enough, J2 is positive. 

The normalized one- and two-body components of 
the vertex functions arc thus 

= - , , V^1 = 0, (54a) 



2myJl+gl^J2 
bY = , , bl' = 0. 



(54b) 



* Neglecting the term with i = 1 brings corrections of relative 
order (m/mi log(mi/m), which tend to zero as the PV 
fermion mass mi increases to infinity. 
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3. Determination of the internal bare coupling constant 

According to our renormalization conditions de- 
tailed in Sec. IIIIBl we calculate the internal BCC 
go2 from the requirement that the w-independent part 
of the two-body vertex function for physical particles 



(i.e. the component fe]""), taken at 



7^, is identi- 



fied with the physical coupling constant g. Since in 
the two-body Fock space truncation 6^^ are constants, 
we immediately get 



.902 



(55) 



It is a well defined, non-pcrturbativc, condition which 
is very convenient to impose in any numerical calcu- 
lation. 

From Eq. ((55|) it follows 



902 



1 - 5^2 



(56) 



The final form of the normalized (and renormalized) 
solution for the vertex function components becomes 



by proper adjustment of the value of fii). Then, we 
express the calculated Fock components through the 
physical coupling constant, normalize the state vec- 
tor, calculate the observables, and after that go over 
to the limit of infinite PV masses in the final result. 
We will demonstrate below how this scheme works in 
practice. 



4- Determination of the external bare coupling constant 

To determine the external BCC one has to consider 
the 3PGF, denoted by G, which is represented, in the 
two-body approximation, by a sum of the two contri- 
butions shown in Fig. 1101 Similarly to the two-body 
vertex function (|36bp . the 3PGF can be decomposed 
in invariant amplitudes: 



S02 




,0 Vl - 9'^J2 ,1 „ 
2m 



g, bl^^O. 



(57a) 
(57b) 



The one- and two-body contributions to the norm of 
the state vector are 



/i = 1 



,9^2, 



.9^2. 



(58) 



As we said above, if fii increases, the quantity J2 
increases too. At some value of fXi we inevitably meet 
the condition J2 = 1, leading to a pole on the r.-h.s. 
of Eq. (inH) , analogous to the well-known Landau pole 
in QED. Further increase of /ii makes gQ2 negative 
and go2 purely imaginary [as well as Eq. (|54ap ]. 
At the same time, the one- and two-body norms Ii and 
I2 become infinitely large, but have opposite signs: Ii 
is negative, while I2 is positive. So one can not get rid 
at this level of the regularization parameters by tak- 
ing the limit fix 00 without formal contradiction. 
A possible way out is as follows. The PV masses play 
an auxiliary role and appear in intermediate calcula- 
tions only, while physical observables must be inde- 
pendent of them (the BCC 502 is not an observable!). 
We do not give any physical interpretation to inter- 
mediate results found with finite PV masses. Below 
we will treat J2 as being a finite quantity satisfying 
the inequahty l — g'^J2 > (formally, this can be done 



FIG. 10: Fermion-boson 3PGF in the two-body approxi- 
mation. 



u{p')Gu{p)^gu{p')[F + B, 



m i/j 

' LU-p 



u{p). (59) 



The invariant functions F and B^^ (scalar form fac- 
tors) depend, under the condition oj-q = 0, on = 
(p' — p)^. Note that the general decomposition ([55)) in 
the Yukawa model is valid in any approximation, since 
the number of independent invariant amplitudes is 
completely determined by particle spins and by sym- 
metries of the interaction. The external BCC is found 
from the requirement 



u[p')Gil[p) = gu{p')u[p) at q^ 



(60) 



Hence, wc must have F{q^ = /i^) = 1 and B^{q^ ~ 
/i^) = 0. In exact calculations, as well as in a given 
order of perturbation theory, we would indeed get 
B^ = because physical amplitudes can not de- 
pend on the LF plane orientation. In approximate 
non-perturbative calculations however it may turn out 
that B^ does not disappear, even at q^ = fi^ . Then, 
according to our renormalization prescriptions, one 
should add to the LF Hamiltonian the counterterm 

lu-a 
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The structure of the counterterm is the same in all ap- 
proximations, while the constant Z3 must depend on 
the Fock sector, in complete analogy with the BCC's. 
We therefore have to calculate both Z31 (for internal 
coupling) and Z31 (for external coupling). In the par- 
ticular case discussed here, i.e for the two-body Fock 
space truncation with the use of the PV rcgulariza- 
tion (with one PV boson and one PV fcrmion), we 
find ~ 0. We therefore get Z31 = Z32 = Z31 = 
^32 = 0. 

We will now concentrate on the calculation of the 
physical scalar form factor F . It can be extracted 
from the 3PGF by the relation 



1 



oF ^ ^ Tr 



(/ + j7i)G(:^ + m) ( !^ - 1 

LO-p 



(61) 



where the variables x and Rj_ are defined by Eqs. ([37]) . 
and r's are given by Eqs. (jC2[) . (jC5|) . and (|C6|) taken 
for M — m. Kinematical relations needed to express 
the scalar products of the four-momenta through x 
and Rj^ can be found in Appendix IC 21 To make the 
integral in Eq. (j63bp convergent, it is enough to reg- 
ularize it by the bosonic PV subtraction only, as we 
did for the calculation of J2 in Eq. (|5ip . Wc can thus 
neglect in the sum the terms with cither i or i' being 
1. According to the solution (j57bp for the two-body 
component of the state vector, T2 = iif = <?. Using 
the solution (|57ap for the one-body component, sub- 
stituting Eqs. ((55a|) and ((63b)) into Eq. ^ and then 
into Eq. ((6T|). wc find 



First, we calculate F. as a function of , for physical 
values < and then find its analytical continuation 
to the non-physical point = IJ-^- 

In the two-body approximation, the 3PGF can be 
represented as 



G = g^Gi + 501G2, 



(62) 



where Gi,2 denote the contributions (amputated from 
the external BCC's) to the full 3PGF from the one- 
and two-body sectors, respectively. Applying the 
CLFD graph technique rules to the diagrams shown 
in Fig. [ini wc have 

= ^(-l)'+'VlV'f [(/ + ^r{+m,0 



(64) 



G2 



X (^+ iJjTi +mi)] 
1 /" ,o „ /"^ dx 



(27r)3 



(63a) 



[p' - fc2j)2 - mj, 



ip- 



(63b) 



where 



^(1) 



= I-.9V2, 



167r- 



(65a) 



(-l)J'a;[2(R^ A)VA2 ~ - xIlj_A + {^x^ - 8x + A)m^ 



c)[i?i 



+ (1 - x)^i'^, + x2m2][(R^ - a;A)2 + {1 - x)^i] + x'^m? 



(65b) 



20 



with A = and -q^. The functions T'^ ^^ , J^^^) 

determine the contributions to the scalar form factor 
F from the one- and two-body Fock sectors. The in- 
tegration over <PR±^ in Eq. (|65b[) can be performed 
analyticaUy and leads to the result 



[4(1 - x)m'^ + 




log 



) 2^x'^m? + (1 - x)n 



(66) 



where 



rj{q^) = y^a;2(4m2-g2)+4(i_^)^2_ 



According to Eq. (j35bp . goi =5, and Eq. ([64)) con- 
tains one unknown parameter, namely, the external 
BCC go2 which is found from the renormalization con- 
dition F(g2 = _ ^ ^Yax^ g^,^ 



.9 = 502(l-5V2)+g.F(2)(^2)^ 
where J2 is defined by Eq. ([55)) and 



(67) 



from various LF diagrams taken into account in our 
calculations. Since 1702 and po2 differ already at order 
of their perturbative expansions, it is enough, to 
clarify the situation, to consider the lowest order per- 
turbative corrections to BCC's. They are indicated in 
Fig. [11] where the outgoing wavy line corresponds to 
an external boson. As it has been pointed out above, 
only the diagrams (a) and (b) are incorporated when 
calculating the internal BCC 502 in the two-body ap- 
proximation. Concerning the external BCC 502, the 
diagrams (c), (d), and (e) are included in addition, 
though they formally correspond to a three-body state 
admixture, with the external field. In other words, our 
calculation of 302 is effectively performed to higher or- 
der approximation, than that of go2- One may expect 
that increasing the number of Fock components in- 
volved into calculations will reduce the difference be- 
tween the internal and external BCC's. 

The diagram (f) in Fig. Ill) containing the fermion- 
antifermion pair intermediate state, does not con- 
tribute, in the two-body approximation, neither to go2 
nor to ^02- For this reason, both BCC's, even being 
expanded in powers of g up to terms of order , differ 
from the BCC found in the second order of perturba- 
tion theory. 

As follows from Eqs. ([64)) and ([67)) . the scalar form 
factor F is given by 



F = I + T^^\q') - 



(2)/ 



(71) 



[4(1 



x)m2 -I- /i2i 



rj(/z2) 



■ arctan 



fix 



The asymptotic value of JF*^2)(^2^ _^ 
^^^HA^^) = -f^log^+5Vi(^) 

OTT^ TO Vto/ 



00 IS 



(68) 



(69) 



where /i is a function of the ratio /x/to, finite at /xi 
00. Substituting Eqs. ([69]) and ([53]) into Eq. ([67]), we 
find for the external BCC: 



1 



902 ^ g- 



log g fi[ — 

OTT^ TO Vto 



167r- 



■log 



Ml 



Vto/ 



(70) 



The internal and external BCC's given by Eqs. ([56]) 
and ([70]) . respectively, differ from each other. We 
have already mentioned in Sec. IIII A 21 that they in- 
deed should not coincide. We can illustrate the ori- 
gin of this difference in a very clear form, by analyz- 
ing contributions to the internal and external BCC's 



It is easy to see that the expression on the r.-h.s. is 
finite in the limit /ii — > 00, though each of the func- 
tions JF(2) diverges logarithmically. Eq. (|7ip coincides 
with the result of the 3PGF renormalization in the 
second order of perturbation theory, in spite of the 
fact that we did not make any expansion in powers of 
the coupling constant. 



B. Application to QED 

1. Determination of the state vector 

The Yukawa model considered in the previous sec- 
tion is a good example of how the general renormal- 
ization scheme developed in this article should be un- 
derstood. We will now apply the method to QED in 
order to address the case of a realistic physical the- 
ory. The perturbative approach to QED was able to 
reproduce multiple experimental data with excellent 
precision. Applying the non-perturbative scheme to 
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two-body Fock space truncation. In this approxima- 
tion, QED and the Yukawa model are very close to 
each other. We repeat below the procedure of find- 
ing the state vector and its renormalization, detailed 
in Sec. IIV A[ for QED in the two-body truncated 
Fock space. Because of similarities between QED and 
the Yukawa model, we will not expose all the steps, 
but concentrate on pointing out the differences and 
demonstrating the main results of our analysis. 

The general structure of the two-body electron- 
photon vertex function in QED has been extensively 
studied in Ref. 01 ■ For the case of the simple two- 
body Fock space truncation, the interaction Hamilto- 
nian in the Feynman gauge is very similar to that in 
the Yukawa model, with the scalar vertex being re- 
placed by the electromagnetic one. Provided the PV 
regularization is used, with one PV boson and one PV 
fermion, the fcrmion self-energy does not depend on 
the LF plane orientation. 

According to Ref. 0, the state vector has the fol- 
lowing structure: 

u{pu)r\u{p) = {m\ - M^)rMPuHp), (72a) 



FIG. 11; Second order perturbative corrections to the in- 
ternal and external BCC's. 



the same object gives us a possibility to test our re- 
sults and to reveal distinctly its main differences, as 
compared to perturbation theory. 

From a purely technical point of view, QED is more 
complicated than the Yukawa model, at least, in the 
following three aspects. Firstly, the structure of the 
interaction Hamiltonian in QED is more involved than 
that in the case of scalar bosons. The considera- 
tion of QED within CLFD, performed in Ref. 0], 
showed that the corresponding Hamiltonian (obtained 
without introducing PV fields) contained specific con- 
tact terms, different from those which appear in the 
Yukawa model. It is not yet clear whether PV fields 
can kill the contact terms in this case. Secondly, 
many-body wave functions in QED have more spin 
components than their counterparts in the Yukawa 
model. Thirdly, intermediate calculations in QED es- 
sentially depend on the gauge condition to constrain 
the electromagnetic field, while the question how the 
Fock space truncation affects gauge invariance is still 
opened. 

The difficulties of treating QED in the framework 
of CLFD, itemized above, are however absent for the 



uiku) Tl^'''ei{k2j) uip) 



= uihi) 



lO-p 



el(k2j)u{p), (72b) 



where e^(fc2j) is the photon polarization vector. The 
scalar functions '^\^ bf , and fej"* differ, generally speak- 
ing, from those entering Eqs. (|36p. The decomposi- 
tion (j72bp with the two matrix components is valid 
only within the two-body Fock space truncation. The 
number of independent components of the two-body 
vertex function in QED depends on the gauge. If the 
Feynman gauge is chosen, there are eight independent 
components in the general case. Under the two-body 
approximation, six of them identically turn into zero. 

The system of equations for the QED vertex func- 
tions is very close to that for the Yukawa model, shown 
in Fig. [9] and given analytically by Eqs. ((38)) and (|39)) . 
Small changes are required, caused by the vector char- 
acter of the photon. Namely, one has to substitute 



502 



6027*^ 



^ 2 



in the Feynman gauge. 



^pi^ip) 



T^2 '^ , and S(p) — > Spj^(p), where. 



Aip') + Bip')^ 
m 



(73) 



The scalar functions A and B here differ from those in 
Eq. ([1^ . They were calculated (for the longitudinal 
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and transversal LFD cutoffs) in Ref. We do not 
need to know their explicit form in the following. 

The technical procedure to solve the system of equa- 
tions for the vertex functions is exactly the same as in 
Sec. IIV AD The solution looks as follows: 



(5to2 = 26^2(2^ - B), 
with A = A{ni'^), B = B{m'^), and 

}hl = 2meo2ai, — 0- 



(74) 



(75a) 
(75b) 



Eqs. (|75ap and (j75bp coincide in form with Eqs. 
and (j46|) . respectively. The constant ai is found from 
the normalization condition. 



2. Normalization of the state vector 

The calculation of the one- and two-body normal- 
ization integrals is analogous to that for the Yukawa 
model. The one-body integral is exactly the same as 
in Eq. (|47p. since it is not sensitive to the boson type. 
The two-body integral is different: 



where 
J2 



h = 4m^aieo2 J2, 



— / (i-R_L Ri^ I dxx 



(76) 



E 



[R\ + (1 -x)^2 ^x2m2]2 



(77) 



The integral J2 at //i — > 00 diverges logarithmically, 
as J2 in the Yukawa model [see Eq. ([55)) ]. but with a 
different coefficient at the logarithm: 



~ I III 1 / TO 9 

•^2 = ^ log T ~ s 

OTT^ TO 47r^ V A 8 



(78) 



where we assigned a finite (small) mass A to the pho- 
ton in order to avoid the infra-red catastrophe. Note 
that J2 is positive in the limit of infinite PV boson 
mass. 

The normalized solution is given by Eqs. (j54p . 
changing 502 by eo2 and J2 by J2. 



3. Determination of the internal hare coupling constant 

Due to the formal coincidence of the state vector 
structures in the Yukawa model and QED, we imme- 
diately obtain from Eq. (|56p : 



p2 



(79) 



The final solution for the vertex function compo- 
nents is 



_ - e2 J, 



2m 



61^ = e 



^ ^/-l = 0, (80a) 
(80b) 



b^i = 0. 



It is in full analogy with the corresponding solution 
for the state vector in the Yukawa model, given by 
Eqs. (Ell). 



4- Determination of the external hare coupling constant 
and calculation of the electromagnetic form factors 

We will establish in this section the relationship 
between the external electromagnetic BCC and the 
physical fcrmion charge and compute the fcrmion elec- 
tromagnetic form factors in the two-body approxima- 
tion. The main difference, as compared to the elec- 
tromagnetic interaction with a system formed by the 
Yukawa interaction, consists in the fact that now the 
photon emitted (or absorbed) by the probe is of the 
same type as photons coming into the state vector. 
Formally speaking, the right diagram in Fig. [10] con- 
tains a three-body intermediate state and must be 
rejected. However, according to our renormalization 
scheme detailed in Sec. IIII A21 one should treat the 
photon connecting the electromagnetic vertex with 
the probe as if it was an external particle, having 
no relation to the contents of the state vector. This 
makes possible the calculation of non-trivial observ- 
ables already within the two-body approximation. 

In CLFD, the spin-1/2 fermion 3PGF (or electro- 
magnetic vertex, i.e. the current matrix element be- 
tween the initial and final fermion states) has the fol- 
lowing general structure [1, [2l| (we have included the 
physical electromagnetic coupling constant e into the 
definition of the vertex): 
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uip')G''u{p) = eu{p') 



2m 



Bi 



A 

uj-p 

r 



pp - 27' 



B2 h i?3 



oj-p 



{Lo-py 



i{p) (81) 



with aP" = - YY)I'2' and P = p' + p. It 

is determined by five form factors, the two physical 
(-^1^2) and the three non-physical {Bi-^) ones. The 
non-physical form factors are coefficients at the spin 
structures which depend on oj. The term proportional 
to Bi is constructed in such a way that it gives zero 
when sandwished between free spinors of momentum 
p' = p. It gives also zero when contracted with ujp. 
Under the condition oj-q = the form factors de- 



pend on 



-q . If the electromagnetic vertex 



is calculated exactly or within a given order of per- 
turbation theory and, moreover, is regularized in a 
rotationally invariant way, the non-physical form fac- 
tors cancel identically, while the two physical form 
factors remain, as it should be. Under approximate 
non-perturbative calculations, the non-physical form 
factors may however survive and plague the electro- 
magnetic vertex by spurious w-dependent contribu- 
tions. The situation here is fully analogous to that for 
the Yukawa model, where the function B,^ in Eq. ([59]) 
is just a non-physical scalar form factor. CLFD allows 
to separate covariantly the physical and non-physical 
parts of the electromagnetic vertex and extract the 
physical form factors from the former. In our case, it 
is enough to contract both sides of Eq. (|8T|) with the 
four- vector ujo- 



u{p' )uj pQP u{p) — eu(p') 



iF2 
2m' 



u{p). 



(82) 

As we see, the three non-physical form factors S1-3 
disappeared, since the contraction of the spin struc- 
tures proportional to them with LUp gives zero. The 
physical form factors can be found by the following 
expressions; 



eFi 



Tr^ji' + m)ujpGP{^ + m)i/j] 



(83a) 



Similarly to Eq. (j62p . we write 



where 



ujpGP = e-mGi + eoiG2, 



(84) 



(85a) 



G2 



(2 



i+i -t-i 



/^2j + TOj')(Z>(^ - /^2i + "^i)r2 



W - k2jY - ml][{p - fc2j)2 - mf] 



(85b) 



Again, to make the integral in Eq. (j85bp convergent, it 
is enough to regularize it by the bosonic PV subtrac- 
tion only. We can thus retain in the sum the term with 
i = i' = 0, neglecting those with either i or i' equal to 
1. According to the solution (|80bp for the two-body 
component of the state vector, T^^^ '^ = ^^'7'' = 67"^. 
Using the solution (|80a[) for the one-body component, 
substituting Eqs. ([85|) into Eq. ([84|) and then into each 
of Eqs. (IMl), we find 



where eoi 



eF2 

= e, 



eo2^}''+eoi.Ff^(Q2), 
eoi^f'(g'), 



(86a) 
(86b) 



eF2 



X Tr 



2(wp)Q2 



i' + m)ujpGP{'^ + m) 



uj-p 



(83b) 



= 1 - eV2, 
J2 is given by Eq. ((77|) . and 



(87) 
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87r3 



(88a) 



- xK^A + (x2 + 2x - 2)m2 - (1 - x)A^] 
[R\ + (1 - a;)At| + a;2m2][(R_L - xA)2 + (1 - a;)At| + x^m^] ' 



47r3 



(88b) 



(-1)^x2(1-2;) 



[i?^ + (1 - X)^^ _^ 2:2^2] [(R_L - xA)2 + (1 - x)/^^ ^ 2:2^2] 



Here /zq = /i = 0. Note that JFp''(0) = e'^J2 which is 
just the norm of the two-body sector. 

The requirement Fi = 1 at = leads to the 
relation 

eFi{Q^ = 0) = e = eo2(l - e^Ja) + eVa, (89) 
so that 



eo2 



(90) 



Hence, the electromagnetic external BCC calculated 
in the two-body approximation is not renormalized. 
Such a coincidence happens not by chance, but reflects 
a general property of the theory. We will discuss this 
fact in more detail below. 

From Eqs (|86p we find the renormalized form fac- 
tors: 



Fi = 

F2 = 



i + .Fp'(g^)-^r'(o), 



(2), 



(91a) 
(91b) 



This result exactly coincides with that found in the 
second order of perturbation theory, though no ex- 
pansions in the coupling constant have been done. A 
similar issue was obtained above for the Yukawa model 
[sec Eq. ((7T|) ]. 

Expressing through the fine structure constant 
a by the relation = Ana and calculating the inte- 
grals we obtain (for simplicity, we decomposed 
the form factors in powers of Q^, up to second order): 



F2 



1 - 



a aQ^ 
27r ^ 127rTO2 



, m 3 
^°^ A - 8 



o(g^). 



+ 0(Q4), (92a) 



(92b) 



The formulas ([921) exactly coincide with the familiar 
perturbative expressions. In particular, F2 at = 
reproduces the well-known Schwinger correction 
to the electron anomalous magnetic moment. As ex- 
pected, the form factors have finite limits, when the 
PV masses tend to infinity. 

The external and internal BCC's calculated in the 
two-body approximation and given, respectively, by 
Eqs. ([501) (|79p do not coincide, as in the Yukawa 
model. This fact has been already explained in the 
end of Sec. IIV A4I In QED however, we encounter 
a new effect, namely, the identity between the exter- 
nal BCC eo2 and the physical electron charge e. It is 
a specific feature of QED, which appears as a conse- 
quence of the Ward identity. To illustrate the situa- 
tion, one may come back to Fig. [TT] Due to the Ward 
identity, the total contribution of the perturbative di- 
agrams (a)-(e) gives zero, and the external BCC eo2 
is not renormalized at all. The same arguments hold 
true for the non-perturbative calculations in two-body 
truncated Fock space. As we will see in the next sec- 
tion, this takes place for truncations of any order. 

One more point is also worth mentioning. Nei- 
ther the perturbative expansion of the internal BCC, 
Eq. dill), at /ii ^ 00, 



-02 



1 



47r m 



coincide with the well- 



nor the external BCC Cg^ 
known purely perturbative formula relating the BCC 
with the physical charge and given in most textbooks 
on quantum field theory: 



1 



a , A -^ 
^log — 



(93) 



where A^ is an invariant cutoff which can be identi- 
fied, in the logarithmic approximation, with /xf . This 
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result is nevertheless quite natural, since Eq. ([M)) . be- 
sides the contributions shown in Fig. [TIT a-e) (which, 
being summed up, give zero due to the Ward iden- 
tity) , takes into account also the effect of vacuum po- 
larization, i.e. the diagram with the electron-positron 
loop on the external photon line, shown in Fig. [TlTf). 
Evidently, our two-body approximation is unable to 
embrace such an effect. The latter can be incorpo- 
rated when Fock sectors containing electron-positron 
pairs arc included into the state vector. 

5. External bare coupling constant. What happens for 
N >2 ? 

We can easily generalize the result to the case 
where an arbitrary number of photons in the state 
vector is considered. We can do that by recurrence. 
In the TV = 3 case (one electron and no more than two 
photons), for instance, the form factors are given by 
the contributions from the diagrams shown in Fig. 1121 
At = we have 




FIG. 12: Electron electromagnetic vertex for the three- 
body Fock space truncation. 



e = eo3^{" + eo2^P'(g' = 0) + eoi^ff '(Q^ ^ o). 

(94) 

The functions ' ' are just the contributions of the 
amplitudes of the three diagrams on the r.-h.s. of the 
graphical equation shown in Fig. [121 For QED, one 
can easily check that Ti^\q^ = 0) = /„, where /„ 
is the contribution of the n-body Fock sector to the 
norm of the state vector. With eoi = eo2 = e we get 

e^eoah+eih+h)- (95) 
Due to the normalization condition, Ii + I2 + I3 = 1- 
We thus have 

e = eo3[l-(/2 + /2)]+e(/2+/3), (96) 



so that 

eo3 = e. (97) 

Repeating the same arguments to the cases of 
N = 4, 5, . . . Fock space truncations, we conclude 
that Eq. ([57)1 holds to all orders in N. Hence, our 
formalism respects the requirement of the Ward 
identity even in truncated Fock space. 



V. APPLICATION TO THE N=3 FOCK 
SPACE TRUNCATION IN A SCALAR MODEL 

We re-examine here our previous calculation po| of 
the state vector in a pure scalar model for the three- 
body Fock space truncation, using the renormaliza- 
tion strategy outlined above. We shall see how this 
strategy should be applied in the context of a non- 
perturbative physical problem where the results can 
not be reproduced in perturbation theory, unlike the 
cases of the two-body approximation in the Yukawa 
model and QED. 

Though all the particles we consider are spinless, we 
will distinguish two types of them, described by the 
free fields ip{x) and x{x) and related to the "main" 
(i.e. analogous to fermions in the Yukawa model or 
QED) and "exchanged" particles. For shortness, we 
will refer to them as scalar "nuclcons" and bosons, 
respectively. The interaction Hamiltonian is 

H'^\x) = -g,^\'~5m\^, (98) 

where the prime at x denotes the fact that it is a sum 
of a physical and a PV component. Since all divergen- 
cies are regularized already by the PV boson, we do 
not need to introduce a PV " nucleon" . The Hamilto- 
nian ((98)) contains the BCC go and the MC 5m'^ which 
must be adjusted to reproduce correctly the physical 
" nucleon" mass m and the physical coupling constant 
9- 

The system of eigenstate equations for the vertex 
functions in the three-body approximation is shown 
graphically in Fig. [T31 Solid and wavy lines corre- 
spond to "nucleons" and bosons. In analytical form 
the system of equations reads 
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SmlTi 
Tl^ - M2 

503ri 



Af2 



(5m2r^(i?^,a;) 



A/2 



+502 X! 



^ (27r)3 

.9o2r^(i?±,a;) 



(l-.T)(si2-A/2) 



dx' 



2x'{l-x-x'){si^3-A'P)'' 



(l-x)(42-A/2) (1-X')(si3-Af2) 



r 



(99a) 



(99b) 
(99c) 



with 



'13 



Rl 



Rl 



(1- 
R'l 



R'\ 



4 



(1 - X') 

(Ri + 



X' 
2 



'123 



(1 - X - x') 



In deriving Eqs. (|99|) wc took into account that fJmf = 
0, since no "nucleon" mass renormahzation occurs for 
the one-body Fock space truncation. As before, the 
index j = relates to the physical boson with the mass 
jjLQ = while j = 1 — to the PV one with the mass 
^i. We also omitted, for simplicity, the superscript 
(3) at each vertex function. Necessary kincmatical 
relations needed to obtain Eqs. (|99p are presented in 
Appendix IC 31 

According to our strategy, the internal BCC 5^02 and 
the MC 5m\ should be determined from the two-body 
approximation. Using the results of Ref. [l^l , we have 

302 = , 'V- . (100a) 
Vl - 9 h 

= -ffo2S(p' = A/2), (100b) 



with 



^2=E 



167r3 



d^Ri 



dx 



x(l-a;)(si2-M2)2 



and 



' 4- 167r3 J ^Jo x(l-x)(s' 



being, respectively, the norm of the two-body sec- 
tor and the " nucleon" self-energy amputated from the 
coupling constant squared. In the limit ^1 —f 00, the 
two-body norm tends to a finite limit, while the self- 
energy diverges logarithmically. 
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FIG. 13: System of equations for the Fock components in 
a three-body pure scalar model. 



The three-body component T3 is expressed through 
the two-body one, r2, and can be eliminated from the 
system of equations ([5^. Since the latter is homoge- 
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neous, we can start with a value of 4>i = Fi/ (m^ — M^) 
equal to one. Dividing Eqs. (|99bp and (I99cp by 503, we 
get a system of inhomogcncous equations for the (non- 
normalized) state vector characterized by the Fock 
components r2.3 = 3/303- Excluding Fi and F3 
from Eqs. ((99|) and going over to the limit M m, 
we obtain a closed inhomogcncous equation involving 
the two-body component only: 



riiR±,x) = 1 



.902 



.902 

167r3 



dx' 
x'{l-x') 



fi{R'^,x') 



(1 



■x'){s' 



13 



™^)(si23 ~"^^) 



where 



R 



+ (1 — x)m^ 



-M7 



(101) 



(102) 



The solution V2 of Eq. (|10ip can be found numer- 
ically for given values of m and g. It is completely 
finite in the limit ^ 00. To get the normalized 
vertex functions, we should calculate the norm iV'^^ 
of the state vector. It can be written schematically as 

Af(3) = AT, + + iVg = 1 + gl^{N2 + 7V3), (103) 

where 1^2,3 are the two- and three-body norms written 
in terms of f 2^3. The normalized components 3 
are thus given by 



1,2,3 



Fi 



2,3 



(104) 



We still have to fix the internal BCC (703 to calcu- 
late the normalized vertex functions. It is found from 
the requirement that the value of the physical (i.e. 
taken for j = 0) two-body vertex function in the 
kinematical point (s* = sj^" = rn^^x*), where x* is 
a given fixed value of x, equals the physical coupling 
constant: 



r^is*,x*)=g. 



(105) 



Let us call Fj the value of Ff in this point. From 
Eqs. pU5)) - (fTU51) we get 



903 



(106) 



which gives 



903 



^[f*Y~g^{N2+N3) 



(107) 



We can see that Eq. (|107p is just a generalization of 
Eq. (jlOOap . since in the absence of the three-body 
component we have Fj = l,iV3 = 0, and N2 = h- 
Using Eqs. (|104p and (|106p . we can represent the nor- 
malized vertex function F|^ simply as 



(108) 



which allows to find it directly from the solution of 
the inhomogeneous equation JlOip . 



without calcu- 



lating explicitly the norm N^^' . The latter is just 
needed to calculate the normalized one-body compo- 
nent <?if = i/Viv(3). 



The substitution of Eq. pU7)) into Eq. pU5)) yields 

1 2 

(109) 



^(3) 



g^N2 + N3) 



Finally, substituting the normalized functions Tf ~ 
(j)f'{m? — NP) and Ff into Eq. (|99ap . we arrive, in the 
limit M TO, at the following expression for the MC: 



5m\ = — 



9 



=.(*) 



d^R 



E 



(27r)3 



dx f{{R^,x) 



2x{l~x) s{2- 



(110) 



The normalized three-body vertex function is found 
from Eq. ()99cp . changing F2 to Ff^. The norms 7V„ 
arc calculated according to 



(27r)3("-i)(n- 1)! 
d^Ri_\_dxi 




(s„ - 



1 



(111) 



where s„ = YM=ii^'i± + "^?)/^' and, in the case of 
the PV regularization, 

f«= E (-i)^^+-+^"-f^-^"-f^-^'-. 
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In such a way, all unknown parameters are expressed 
through the solution of the inhomogcneous equa- 
tion (fTOTI) . 

The application of this strategy to the three-body 
Fock space truncation in the Yukawa model (one 
fermion and two scalar bosons) is in progress. 



VI. CONCLUDING REMARKS AND 
PERSPECTIVES 

We have presented a systematic strategy to calcu- 
late physical observables in CLFD, when Fock space 
is truncated. This requires to implement an appro- 
priate renormalization scheme in a non-perturbative 
framework. Within CLFD, we have shown how to fix 
the counterterms and the BCC's of the elementary 
Hamiltonian in a consistent way. As a check of our 
formalism, we treated the Yukawa model and QED in 
the two-body approximation, as well as a pure scalar 
model in the three-body approximation. We were able 
to recover, for the first time, the standard renormal- 
ization of the electromagnetic charge according to the 
Ward identity, without any perturbative expansion. 
This shows that no divergences are left uncancelled 
in the above calculations. Moreover, the first correc- 
tion to the electron anomalous magnetic moment (the 
Schwinger correction) is recovered analytically. 

Our results have been made possible because of the 
following three important features of the formalism: 

(i) First we can extract the physical part of the 
two-body vertex function in the point s = m^, 
in order to identify it with the physical coupling 
constant. This part is explicit in our formalism, since 
it should be independent of the LF plane orientation 
determined by the four- vector lo. 

(ii) The counterterms and the BCC do depend on 
the Fock sector, in order to cancel all divergences. We 
give a well-defined systematic procedure to calculate 
them. 

(iii) We have to distinguish two types of BCC's: 
those of the first type are used to calculate the state 
vector itself, by means of the eigenstate equation, 
while the BCC's of the second type are responsible 
for describing interactions of the constituents with 
external particles (electromagnetic probes, for exam- 
ple). As we increase the number of Fock sectors, the 
two BCC's of different types are expected to converge 



to the same limit (for finite PV masses). 

Our results are very encouraging in the perspective 
of doing true non-perturbative calculations of bound 
state systems in a field theoretical framework. In case 
of success, they may become a real alternative for lat- 
tice calculations. 
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APPENDIX A: LIGHT-FRONT 
HAMILTONIAN IN THE YUKAWA MODEL 
WITH THE PAULI-VILLARS 
REGULARIZATION 

We shall find in this section the interaction LF 
Hamiltonian incorporating PV fields. The method- 
ology is the same as in Ref. 0]. For this reason, we 
will not give here the detailed algebra with all explicit 
intermediate results, but outline only the main steps 
of the procedure. 

First, substituting the Lagrangian pU|) into Eq. ^ 
and keeping in mind that now Yi runs through the 
whole set of the physical and PV fields, we construct 
the corresponding energy-momentum tensor and then 
find the four-momentum operator in terms of the 
Hcisenbcrg fields. In order to have the interaction part 
of the four-momentum operator in Schrodinger repre- 
sentation, we should find the constraints which con- 
nect different field components at the same time and 
distinguish the independent components. We then ex- 
press the four-momentum operator through the inde- 
pendent field components only and change after that 
all the Heisenberg operators by the free ones. 

It is convenient to perform calculations in the ref- 
erence frame where the four-vector co'' has the com- 
ponents (1,0,0,-1) [the four-vector is (1,0,0,1)].^ 



^ It is not necessary to fix any particular reference frame. All 
subsequent algebraic manipulations, in principle, can be done 
in explicitly covariant notations, however, the calculations 
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Under this condition the LF " time" is = + , 
while the "coordinates" are x~ — x'^ — x^ and x-*- = 
{x^ , x'^). Analogously, we introduce the plus-, minus-, 
and transverse components for any four- vector. 
The equations of motion for the fcrmionic fields are 



(i-f^dr^-m)^ = -(.90$' + (5™)*', (Ala) 
{iYd^-mi)'i>pv = (go*' + (5771)*', (Alb) 

Splitting the bispinors into the so-called plus- and 
minus-components (they should not be confused with 
the plus- and minus-components of four vectors) 



(A2) 



Pp through them. After that, finding this operator in 
Schrodinger or interaction representation reduces to 
simple changing the Heisenberg fields ^pvi ^^ 

and ^pv by their free counterparts 'ip^^\ ^pyi f-nd 
ifpv- Splitting, as in Eq. ([5]), the four-momentum op- 
erator into the free and interaction parts, and keeping 
in mind that ui^ = {luq + uJ3)/2 — 1, u}± ~ —u)-^ ~ 0, 
w_ = {luq — uj3)/2 = 0, we obtain from Eq. 



PI 



Sx^ 



dx 



ipv 



(A6) 



Calculations performed according to the rules given 
above yield the formula pT|) for the interaction Hamil- 
tonian. 



where 



M/W=A(±)vi/, A(±) 



:(1±7V) 



(A3) 



APPENDIX B: CLFD GRAPH TECHNIQUE 
RULES 



(and similarly for 5'pv), we can express, by means of 
Eqs. (|A1|) . the minus-components through the plus- 
ones: 



J L 



-,/(+) 



{-) 

PV 



-(.90$' + (Sm)\E' 



/(+) 



7 

id+ 
+ (5o*' 



(+) 

PV 



SmW 



(A4a) 



(A4b) 



The operator l/(ii9+) acting on an arbitrary function 
/ of the coordinate x involves integration over the 
minus-component of its argument: 



The graph technique rules for CLFD have already 
been given in Ref. Q for the calculation of the 5*- 
matrix in terms of the operator —Hint- Since the 
eigenstate equation ([5^ is also expressed in terms of 
the same quantity, these rules are applicable in our 
case too. 

Because of the particular conservation law ((29)) , it is 
convenient to represent, for the calculation of CLFD 
diagrams, the momentum toTn by a separate line, in 
addition to the ordinary particle lines. This is the 
so-called spurion line |6[. For practical calculations, 
one may thus represent the vertex function, as well as 
any elementary vertex, by the diagrams indicated in 
Fig. 1141 The spurious are shown by the dashed lines. 
The momentum conservation laws associated to the 



dy e{x -y )f(y ), (A5) 



where e is the sign function. 

Eqs. (jA4p are constraints, since they do not include 
the "time" derivatives d~ and connect field compo- 
nents at any time. We see that the minus-components 
of the spinor fields are not independent in the sense 
that they are expressed through \E'(+\ ^pv' ^^'^ 
^pv- Taking the four latter quantities as independent 
fields, we can express the four- momentum operator 




(a) 




FIG. 14: Convenient representation of the vertex function 
(a) and a typical elementary vertex (b) in terms of the 
spurion line. 



would be more lengthy. 
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vertices (a) and (b) in Fig. [T3] are, respectively 
and 

S^'^^p + q-p' +u{t' -t)]. 

We emphasize that spurions are not considered as true 
particles. They serve simply for a convenient repre- 
sentation of the departure of intermediate particles 
from the energy shell. 

Appearance of PV fields changes the CLFD graph 
technique rules, as compared to those listed in Ref. [0]. 
However, one does not need to re-derive the rules 
from the very beginning in order to incorporate the 
changes. Indeed, PV particles can be considered as 
new fermions and bosons having their own masses and 
negative norms. Hence, they participate in the con- 
struction of CLFD diagrams on equal grounds with 
physical particles, simply increasing the number of di- 
agrams. The analytical rules to calculate amplitudes 
of the latter remain almost the same, excepting the 
expressions for particle propagators. Since physical 
and PV fields represent two different groups of parti- 
cles, they can not pair between themselves, i.e. any 
line in the diagrams relates to either a physical or a 
PV particle. Pairings within each group are allowed. 
PV propagators correspond to the PV particle masses 
and also differ by a sign from the physical ones. In 
momentum space the propagators are 



-iye{uj-p)5{p^ 



for bosons and 



(-l)X;^ + 7n,)e(^-p)'5(p'-m,2) 

for fermions. The indices i and j describe the particle 
type: i = and j = correspond to physical particles, 
while i = 1 and j = \ relate to PV ones. 



APPENDIX C: KINEMATIC AL RELATIONS 

1. One- and two-body kinematics 

To calculate the amplitudes of the CLFD diagrams 
contributing to Eqs. ((38)) and p9)) . we need to express 
the intermediate momenta and their contractions with 
the external momenta through the variables , R'l , 
X. and x' . 



The momentum conservation laws in the vertices 
lead to the following equalities: 



pu = P + t^Tl, 

fcli = P + - k2j, 

Ki' = P + ^T2- k'2j- ■ 



(Cla) 
(Clb) 
(Clc) 
(Cld) 



Squaring the first two of these equations and taking 
into account that p^ ~ AP, p\^ ~ mf, p'\^, — mf,, and 



0, we get 



2lo-p 



2uj-p 



(C2) 



Since, by the definition (|37|1 . LU-{k2j —xp) ~ 0, we have 



{k2j - xpY 



-(k: 



. This allows 



to represent the scalar product p-k2j as 



p-h 



■K2j 



Rl + M^x^ 
2x 



(C3) 



For the scalar product p-k'^j, we analogously find 



p-k'2 



2j' 



2x' 



(C4) 



Squaring Eq. (|Clc[) under the condition k\^ = m^, we 
obtain with the help of Eq. (jC3|) : 

mf ~ {p- 



T2 



k2,r 



2uj-{p— k2j) 



2u}-p 
Analogously, 



1 



2uj-p 



1-x 



1 - x' 



(C5) 



(C6) 



On the other hand, since {p + cur.^ 



s^2 I we have 



* 7 



+ 2{lu-p)t^ 



R'l + ml 
1-x' 



(C7) 



The above equalities are enough to calculate the self- 
energy (j4T|) and to obtain the system of equations (|43|) 
for the scalar components of the vertex functions from 
the initial system of equations 
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2. Two-body contribution to the three-point 
Green's function 

To calculate the two-body contribution to the 
scalar and electromagnetic form factors from the ver- 
tices (j63bp and (|85b|) , we need to express the contrac- 
tions of the four- vector k2j with the external momenta 
p and p' through the integration variables Rx and x. 
The contraction p-fc2j is given by Eq. (|C3p . In order to 
find the scalar product p'-k2j, we introduce the four- 
vector R' = k2j — xp' and consider its properties. As 
Lo-p ~ i^-p' , we have uj-R' = and. hence, 



r" 



~R'l 



-(k2j_L - xp± - xq±) 



= (Ri - xAf. (C8) 

On the other hand, 
R'^ = {k2j - xp'f = pi] - 2x{p'-k2j) + x^m^. (C9) 



Comparing the r.-h.s.'s of Eqs. (jC8[) and (|C9[) . we fi- 
nally get 

(R^ - xAY + + x'^m^ 
p'-k2, = . (CIO) 



3. Three-body kinematics 

To cast the equations for the vertex functions in the 
three-body approximation to the form (|99p . wc should 
define due kincmatical variables. For this purpose, we 
consider a set of three four- vectors 



The substitution of fc„ = Rn + XnP into Eq. (jC12p 
gives 



and, hence, 



(CM) 



(C15) 



Because of Eqs. (jC13|) and (|C15[) . only two of the 
three R_l's and two of the three x's are independent. 
We will choose as independent variables those for the 
bosons, while the "nuclcon" variables are expressed 
through them as 



Rii = -(R2± + R3±), xi = l-x2-x^. (C16) 

We also need to know the contractions of the four- 
momenta p and fc„ among themselves. Since for all n 
we have uj-Rn = 0, i?^ = —R^^, then, in full analogy 
with the derivation of Eq. (|C3p . 



p-k„ 



(C17) 



where m„ is the mass related to the particle n. From 
the condition uj-Rn = follows 



(fcrii ^ XniP)'{kn2 ~ Xn^p), (CIS) 



which finally yields 



(Cll) 



where n = 1, 2, 3 and Xn ~ {i^-kn) / [oj-p). The 
four-vectors fc„ are just the four-momenta of the con- 
stituents. For shortness, we do not pay attention to 
the particle type (physical or PV one). The momen- 
tum conservation law reads 



(C12) 



where ujt^ is the spurion momentum in the three-body 
state. Contracting both sides of Eq. (jC12p with w, we 
have 



= 1. 



(C13) 



h -h 



(C19) 



-|--^(i?^^j_ -f- ml j - R„j_lR„2^. 

It is also convenient to define the invariant energy 
squared of the three-body state: 



Sl23 



(C20) 



Omitting the subscript 2 at R2± and X2, and chang- 
ing Ra^ R'l, x' everywhere in the above 
formulas, we transform the equations for the vertex 
functions, obtained from the diagrams of Fig. [121 to 
those given by Eqs. ([M]) . 
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